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I. INTRODUCTION 

The limit of the large number of colors Nc oo and the 1/Nc expansion in QCD are powerful methods of theoretical 
analysis which provide us not only with a qualitative understanding of nonperturbative phenomena of QCD but also 
sometimes lead to interesting quantitative results. The history of the applications of the 1/Nc expansion in QCD is 
almost as old Q, |^ H, Q as QCD itself and the roots of this method in statistical and many-body physics are even 
deeper. 

In spite of this long history some quantities important for the phenomenological applications of QCD have for a 
long time remained out of the scope of the large- iVc community. In this paper we construct the 1/Nc expansion for 
the baryon wave function. The light-cone version of the baryon wave function (distribution amplitude) plays a crucial 
role in the QCD analysis of hard exclusive phenomena [EISSQ- 

Apart from the importance for practical applications, the baryon wave function is a quantity with a very interesting 
structure of the 1/Nc expansion. The traditional methods of the 1/Nc expansion usually deal with quantities that 
have a power lavge-Nc behavior. For a vast amount of baryon properties (masses, mass splittings, magnetic moments, 
form factors, structure functions, etc.), the large-A'^c behavior can be described by a power counting of Nc. On the 
contrary, the large- A'c analysis of the baryon wave function requires the development of methods that can be applied 
to quantities which depend on Nc exponentially. 

The main part of this paper is devoted to the light-cone case which is especially interesting for practical applications. 
However, for simplicity we prefer to start the analysis of the large-TVc behavior from the "covariant quark wave 
function" of a baryon \B), 

*f/l^l)(/2S2)...(/iV,Sivj(^l'^2,- •■,ZArJ EE .^£^^_^^JO\T [qc^nsAzi)qC2f2S2{z2) ■•■ (^wj] \B) . (1.1) 

Here qcfs{z) is the quark field taken at the space-time point z, and c, /, s are color, flavor and spin indices, respectively. 
The Levi-Civita tensor £c-i...cn^ is totally antisymmetric in color indices. We use the notation 4" with a hat in order to 
distinguish the covariant wave function 1)1. l|l from the light-cone case, which is the main subject of the paper. Because 
of the time ordering T and Fermi statistics of quarks, the function ^ff-^^sx) (/« j^'^i' ■ ■ ■ '-^^'c) symmetric with 
respect to permutations of {zk, fk, Sk)- Special care should be taken about the local gauge invariance. We assume 
an insertion of Wilson lines on the RHS of Eq. (|l.l|l in such a way that function \£'^ becomes gauge invariant but 
still remains symmetric with respect to permutations of {zk, fk, Sfe) (for example, we can choose straight Wilson lines 
connecting each quark field with the "center of mass" of points zi, Z2, . . . , zn^ and antisymmetrize the Nc color indices 
at this center point). 

Our aim is to understand the large- TVc behavior of the baryon wave function (|1.1|) . But if we try to apply the 
1/A^c expansion to the function (|l.ll) directly, then we immediately meet a problem. This function depends on Nc 
arguments. From the mathematical point of view it hardly makes sense to speak about the large-TVc asymptotic 
behavior of a function depending on Nc variables. For a consistent treatment of the large- A^^ limit we need some 
intermediate object which 
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1) contains the same information as the wave function w , \(zi, zo, 

2) allows to take the large- -/Vc limit and to construct the 1/iVc expansion. 
To this aim we introduce the generating functional 



^b(5) = XI / / dz2---J '^^We5/isi(2l)ff/2^2(^2)--.g/„,s„J^ArJ^['(/isi)(/2S2)...(/«,s«j(^l'^2^ (1-2) 

fkSk 

depending on an arbitrary "source" function gfs{z). Since the function is totally symmetric in permutations of 
{zk, fk,Sk), at finite Nc we can restore the wave function if we know the functional ^^(ff): 



B 

(/lSl)(/2S2)---(/jVcSNc) 



1 



{zi,Z2,...,ZnJ = -TTT 



n 



(1.3) 



What can we say about the large- TVc asymptotic behavior of the functional $5(9)? In principle, it is a dynamic 
problem and without solving large- A'c QCD we cannot give an absolutely reliable answer. Nevertheless on general 
grounds it is natural to expect the behavior of the type 



,{g) = {k-Ab{9) exp [N,W{g)] } [l + 0{N-')] 



(1.4) 



Here the exponential behavior is controlled by the functional W{g). The pre-exponential factor Asig) may be 
accompanied by a power term N^^ . In this paper we present a number of arguments in favor of the large- A^c behavior 
()1.4|l . We shall also establish some interesting properties of functionals W{g) and Asig)- In particular, we shall show 
that the functional W{g) is universal for all low- lying baryons (see Sec. IIIII for details). 



II. MAIN RESULTS 
A. Notation 

In this paper we deal with the l/Nc expansion for similar but different versions of the baryon wave function. In all 
cases we introduce generating functionals analogous to ^3(9) lll.2|l and describe their large- A^c asymptotic behavior 
in terms of functionals similar to W{g) and Asig) in Eq. (|1.4|) . In order to avoid confusion we use a slightly different 
notation for various cases: 

• covariant baryon wave function (|1.1|) and the associated functionals $s(g), W{g), Asig), 

• baryon wave function V'^ (|4.2|l in a toy quark model and functionals (psig), w{g), asig), 

• baryon light-cone wave function (distribution amplitude) (|5.5|) and functionals ^3(9), W{g), Asig), 

• asymptotic baryon distribution amplitude and other wave functions diagonalizing the evolution kernel (|9.5|l . 
and functionals $B(g), Wo{g), Asig)- 

On the other hand, the general structure of the 1/Nc expansion in all these cases is almost identical. Therefore in 
the rest of this section describing the main results of the paper, we use the same notation ^3(9), W^(ff), ^3(9) for 
all above cases. The precise meaning of these functionals will be obvious from the context. 



B. General structure of the 1/Nc expansion and universality of the functional W{g) 

One of the aims of this paper is to check the consistency of the large- A^'c structure (|1.4|) of the functional $5(5)- 
The traditional machinery of the large- A'^c analysis includes several approaches: 

• equivalence of the large-A'c coun ting in QCD and in various models imitating QCD 0, SH HJl E3, 0| 
(nonrelativistic quark model [13 . Il3l IT^. Skyrme model 0, 0|, chiral quark-soliton model |lS Il9l l20|. 
etc.), 

• direct counting of Nc orders in perturbative Feynman diagrams Q, |M S 1 
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• methods based on the spin-flavor symmetry and consistency condition |^ 0, |^ . 

In Sec. IIVI we compute the functional $B(.g) in a toy quark model and show explicitly that the large- iVc behavior 
of this functional has the structure H1.4|l . Using this toy model, we describe the saddle point method which will be 
later applied to the analysis of the large- A'c behavior of the baryon distribution amplitude in QCD. 

Now let us turn to the role of perturbative QCD in the analysis of the large- iV^ behavior of the baryon wave 
function. Perturbative Feynman diagrams allow us to estimate the Nc order of various nonperturbative quantities 
even if the perturbation theory does not work for these quantities. Unfortunately this method cannot be applied to 
$5(5) directly. Indeed the functional ^3(9) depends on Nc exponentially so that arbitrarily high orders of Nc can be 
met in perturbative Feynman diagrams. Nevertheless perturbative QCD provides a very important tool which allows 
us to check the large- A^c structure of $b {g) described by Eq. I|1.4|l — we mean the evolution equation for the baryon 
distribution amplitude. The main part of this paper is devoted to the evolution of the baryon distribution amplitudes 
at large Nc- 

We stress that in this paper we are not interested in the physical applications of the evolution of the baryon 
distribution amplitude. The 0{Nc) growth of the nucleon mass at large Nc changes the physical picture of hard 
processes in large- A^c QCD. In fact, in practical applications there is no need in the large- A'c approximation for the 
evolution. 

Instead, we would like to use the evolution equation as a consistency test for the large- iVc structure p.4|l . In some 
sense the evolution equation plays the same role in the analysis of the exponential large- iVc behavior as the perturbative 
Feynman diagrams in the analysis of quantities with the power large- iVc asymptotic behavior. In physical applications 
to hard processes the evolution equations are used for the exponentiation of large logarithms of the hard momentum 
log Q. In the large- TVc analysis of the baryon distribution amplitude, the evolution equation is helpful due to its ability 
to exponentiate Nc- 

An important property of the l/Nc expansion is the universality of the functional W{g): this functional is the same 
for all low-lying baryons [including 0{Nc^) and 0{Nc) excited resonances]. Moreover, the same functional W{g) can 
be used for the quark wave functions of baryon-meson scattering states. This universality property is described in 
detail in Sec. IIIII All methods used for the analysis of the large- A'c behavior in this paper (toy quark model, evolution 
equations, asymptotic limit) confirm the universality of the functional W{g). 

In contrast to the W{g), the functional Asig) appearing as a pre-exponential factor in Eq. 11.4|l depends on the 
baryon (baryon-meson state) B. We find interesting factorization properties of the functional Asig) which are briefly 
described in Sections lIII DI andlXl The complete analysis of these factorization properties would require the techniques 
based on the large- A^c contracted SU {2Np) spin-flavor symmetry [3, flOl I21L 122. |23| . Here we do not touch upon this 
method, leaving it for a separate paper |2J| . 



C. Nonlinear evolution equation and its Hamilton-Jacobi structure 

In Sec. IV Al we introduce the light-cone version of the functional $3(3) describing the baryon distribution amplitude 
and define the corresponding functionals W{g) and Asig) using the large-A''c representation H1.4|) . In Sec. IVII we 
derive the evolution equations for W{g) and Asig)- The resulting evolution equation (|6.13l) for the functional W{g) 
is nonlinear. This nonlinear evolution equation has the Hamilton-Jacobi form (Sec. IVID|I . The asymptotic limit 
of the large normalization scale /i ^ cxd is studied in Sections IVIII IVIIII In Sec. IIXI we solve the problem of the 
diagonalization of the evolution of the baryon distribution amplitude at large Nc and compute the corresponding 
anomalous dimensions in two orders of the 1/A^c expansion. 



D. Conformal symmetry and integrability 

An essential part of this paper is devoted to the evolution equation and to the problem of the diagonalization 
of anomalous dimensions of baryon operators at large Nc- At Nc — 3 this problem has attracted much attention 
[25I I28L I29I Isof and was studied with various methods ranging from straightforward numerical attacks to the 

discovery of the integrability of the evolution equation for baryons with helicity 3/2 31. .32j. 

The conformal symmetry of the evolution equation is known to play an important role in this problem |^ I27I IsM 
[si*, '35', 'sg'I . It is well known that the large- A'c limit is often helpful for solving problems inspired by high-energy QCD 
[37, 38] where the dynamics is determined by the gluon sector. In these problems the large- A^c limit leads to the "spin 
chains" with the interaction only between the next neig hbors |39|. The case of the quark distribution amplitudes of 
baryons is different. In this case we have equally important pair interactions between all Nc quarks. We derive a 
variational equation for the asymptotic limit of the functional W{g) and solve this equation analytically in Sec. IVIIII 
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Technically the solution is found in terms of a successfully guessed ansatz but the role of the conformal symmetry 
standing behind this ansatz is rather visible. Another problem solved analytically in this paper is the diagonalization 
of the anomalous dimensions of baryon operators of leading twist in two orders of the \/Nc expansion \0{Nc) and 
0{N^)]. The calculation of these anomalous dimensions in the next-to-leading 0{N^) order is based on the solution 
of rather nontrivial functional equations. The analytical solutions of these equations are found in Sec. IIX CI 



E. Phenomenological applications 

On the phenomenological side one is certainly interested in the calculation of the functionals W{g) and Asig) 
appearing in the 1/Nc expansion H1.4|l . This challenging problem belongs to the same class of difficulty as solving 
large- iVc QCD. Although we do not know the functionals W{g) and ^^(g), in this paper we establish several properties 
which may be of interest for the phenomenological applications. One of them is the universality of the functional 
W{g) which allows a unified description of different baryons and meson-baryon scattering states. 

The large- iVf, arguments are sometimes used j40L l4ll | for the justification of various approximations in phenomeno- 
logical models of baryons. Our results impose certain theoretical restrictions on these models. As was explained above, 
it is probably impossible to construct a systematic expansion for the baryon wave function directly, without 

using auxiliary objects like the generating functional (|1.2|) . Therefore the models compatible with the l/N^ expansion 
can access only the weighted integrals of the baryon wave function rather than the wave function itself. One should 
keep in mind that the one-to-one correspondence H1.2|l . between the wave function and the functional ^sig) 
holds only at finite N^- We cannot compute the wave function at finite Nc by applying the variational derivatives 
of Eq. H1.3|l to the large- iVc asymptotic expression ()1.4|1 for $5(3). 

Another theoretical restriction imposed on the phenomenological analysis by the systematic 1/Nc expansion is the 
necessity to work with the A^c-point baryon wave function [via the generating functional ^^(g)]. The arguments of 
the 1 /Nc expansion cannot be used in models explicitly dealing with baryons "made of three quarks" . Although one 
may think that working with the 3-point baryon wave function is equivalent to a partial resummation of higher 1 /N^ 
corrections, the restriction to the 3-point baryon wave functions is incompatible with the systematic 1/A^c-expansion. 

Obviously the universality of the exponential large- A^c behavior for all low-lying baryons is not sufficient for serious 
practical applications. One also needs a good control of the pre-exponential non- universal functionals Asig)- The 
factorizable structure of these functionals is studied in detail in Ref. [2J| where an interesting relation connecting the 
cases of nucleon and A resonance is derived. 



III. UNIVERSALITY OF THE GENERATING FUNCTIONAL W(g) 



A. Universality of W{g) for baryons 



As was already mentioned, the functional W{g) appearing in the large-A^c decomposition (|1.4|l is universal for all 
low-lying baryons. This property is rather general: it also holds in the case of the light-cone distribution amplitude 
which is the main subject of the paper and in the exactly solvable toy quark model that will be considered in Sec. II VI 
In order to be specific, we describe this universality property for the case of the covariant baryon wave function 
keeping in mind that all statements made in this section are also valid in the light-cone case. 

First we must specify the baryons for which the universality of the functional W{g) is expected. According to the 
standard picture of baryons in large- A^c QCD with Nf — 2 flavors, the lowest baryons have equal spin J and isospin 
T Ellll: 



T = J 



1 3 

2 ' 2 ' • ■ 

0,1,2,. 



for odd Nc , 
for even A^^ . 



(3.1) 



The masses of these baryons have a. 1/Nc suppressed splitting: 



MT=.j^diNc + da + Nc ^ 



Sl\+d>^[j{J+l) +0(A^-") 



7(2 



(3.2) 



These baryons belong to the same representation of the spin-flavor symmetry group |^ ITol I211 123 . |2^ which becomes 
asymptotically exact at large Nc- As a consequence, the functional W{g) introduced in Eq. (|1.4|l is the same for all 
low- lying baryons (|3.1|l . The dependence on the type of baryons B appears only in the pre-exponential factor AB{g)- 
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Moreover, at large Nc one can also consider higher baryon resonances with the excitation energy O(N^). Their 
mass spectrum is described by the formula 

M = diN, + (^do + Y, ^^"^ j + OiN^^) , (3.3) 

where = 0,1,2, .. . are integer numbers, and the cui are N^, independent coefficients [we have omitted the higher 
0{N~^) correction depending on additional quantum numbers]. At large Nc » rii the functional ^{g) for these baryons 
also has the structure H1.4|l with the same "universal" functional W{g) but with baryon-dependent pre-exponential 
fimctionals Asig)- 



B. Baryon-meson scattering states 



Strictly speaking, the existence of the states (|3.3|l depends on the mass thresholds which control the decays of excited 
baryons into lower baryons and mesons. If these decays are possible, then the baryons H3.3|l have a parametrically 
large width 0{N^). In this case instead of large- width baryons it is better to work directly with the scattering states 
containing one baryon and one (or several) meson jij, 13 E3. 1^ ■ We still can use expression (|3.3(l for the 

description of the energy of these scattering states identifying parameters cui with the energies of single mesons in the 
rest frame of the 0{Nc) heavy baryon. 

For these baryon-meson scattering states one can also define the wave function (|1.2|) with the large- iVc behavior 
(|1.4(l . This representation H1.4|l will contain the same "universal" functional W{g) as before (we assume that the 
number of mesons in the scattering state is kept fixed in the limit of large Nc). 

The spectrum of applications of the universality of the functional W{g) can be extended. For example, instead 
of the scattering states with one baryon and on-shell mesons we can consider the matrix elements with additional 
insertions of off-shell quark-antiquark color-singlet currents 



Jk{y) = X! ic{y)Skqc{y) 



(3.4) 



where Sk is some spin- flavor matrix. If the number n of these insertions is kept flxed in the large- A'c limit, then the 
large- A'c behavior 

-^eci...c„, XI / / dz2---J dzN,gf,sAzi)9f2S2iz2)---gh^sr,AzMj 

/fcSfc 



x(o|r. 



.k=l 



(lcihs^{zi)qc2hsAz2) ■ ■ ■qcMjN^sr.S^N^ \ \B) 



N? AB,{j,(^y,)}{g)e^V NcW{g) 



(3.5) 



is described by Eq. (|1.4|l with the same functional W{g) as in the case of the baryon wave function. 

Let us consider QCD with the exact SU{2) flavor invariance. At flnite Nc the functional '^sig) is invariant if we 
simultaneously perform a flavor rotation of the quark flelds and of the baryon state. But due to the universality of 
the functional W{g) the flavor rotation of the baryon state is not needed. Thus we conclude that functional W{g) is 
invariant under flavor rotations of g. Similarly, the universality of W{g) for baryons with different polarizations leads 
to the invariance of W{g) with respect to the spacial rotations of g [including both the rotation of the space argument 
z and the transformation of the spin index s of gfs{z'^, z)]. 



C. Universality and SU{Nf) symmetry 

Note that the universality of W{g) is based on certain assumptions about the large-iVc behavior of the quantum 
numbers of baryons. For example, in the case of the SU{Nf) flavor symmetry with Nf — 2, the functional W{g) is 
universal only for those baryons whose spin J and isospin T do not grow with Nc'. 



I,J = 0{N°J. 



(3.6) 
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In the case of Nf = 3 we must also keep the strangeness S bounded in the Hmit iVc — > oo 

S = 0(7VO) , (3.7) 

if we want to have the universahty of W{g). 

D. Factorization of the pre-exponential factors 

In contrast to the universahty of the W{g), the functional ^s,{j^(j,^)}(.g) depends both on the type of the baryon and 

on the currents JkiUk)- An important property of the functional ^B,{jfc(j,fc)}(5) is its factorization. In an oversimplified 
form this factorization is 

,{Jk(yk)} i9)=Aoi9)Y[[^^i9)Y'Yl^VJA9)] , (3.8) 

i k 

where the functional Ao(g) corresponds to the wave function of the lowest baryon, the factors £,i{g) are associated 
with the elementary uji excitations in Eq. (|3.3|l . and the factors fij^ig) correspond to the insertions of currents Jk in 
Eq. H3.5|l . The precise expression for ^s,{j^(j,^)}(.g) must also contain functional factors which come from the zero 
modes corresponding to the spin-flavor rotations. The role of zero modes is discussed in Sec. IX CI and in Ref. [2^ . 



IV. TOY QUARK MODEL 



A. Model 



We want to compute the analog 0b (.g) of the functional $_b(5) in a simple toy model and to show explicitly 

how the large- A^c asymptotic behavior (|1.4|l appears. This can be done in two ways. First we compute the functional 
(f>B{g) exactly and take the large- A^'c limit of this exact result. In Sec. IIV CI we describe another approach based on 
the saddle point method. The results obtained in the toy model will also be used in Sec. IVIII CI where we study the 
asymptotic solutions of the nonlinear evolution equation. 

Our toy quark model deals with baryons made of Nc quarks. We assume that quarks with flavor /, spin s and color 
c are pinned at one point so that the spacial motion is ignored. We consider the case of Nf = 2 flavors. 

In this simple model, the spin J of a baryon coincides with its isospin T 

ri,|,...,^ foroddiV,, , . 

^ ^ " 1 0,1,2,..., for evenTV,, 

and the quark spin-flavor wave function of a baryon can be written in the form 

fe=i 

Here D-'^^,{R) are Wigner functions for the group SU{2). The integral runs over the SU{2) matrices R with the 
Haar measure normalized by the condition 

' dR^l. (4.3) 
The normalization coefficient ctn^ is chosen in Eq. H4.2|l so that 

I 2 

X! i'^r/f'fK/2'*2)...(/jves«j (4-4) 

fkSk 

This normalization constant is computed in Appendix IXI [see Eq. (|A16p ] 
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B. Generating function (f'TTsJ^fig) 

By analogy with Eq. H1.2|l we introduce the generating function 

/fcSfc 

= CTN^ J dRDl^,^{R-^) [Tr {Rg'^)]''' , 

where gfs is an arbitrary matrix. The superscript tr in g*'' stands for the matrix transposition. 
In Appendix 1X1 we compute (l){g) at finite Nc [see Eq. (|A25ll ] 



(4.6) 



CTNa 



{detg) 



1/2 



(4.7) 



Here g (det 5)"^^^ is an SL{2, C) matrix. Therefore the Wigner function D^^ should be understood in the sense of 
the complexification of SU{2) to SL{2,C). 
At large Nc we find from Eq. I|4.5|l 



The large- iVc asymptotic behavior of (f'TT^Jsig) has the standard form (|1.4|l 

^TT.M = {K^aTTsjAg) exp [Ncw{g)]} [1 + 0{N-')] 

with 



wig) = -ln(2det5) 



(4.8) 
(4.9) 

(4.10) 



i^B = -3/4. 



(4.11) 



aTTaJsig) 



1/4 



\/2r + id: 



T3J3 



(detf/) 



1/2 



(4.12) 



We see from Eq. (|4.1U|) that function w{g) is independent of T = J, T3, J3. This is a manifestation of the "universality" 
property of w(g) which was discussed in Sec. IIIII 



C. Saddle point method 

Although the above results H4.1()|l and H4.12(l for w{g) and arTsJsig) can be read directly from the exact finite- A"c 
expression H4.7|l . it is instructive to derive these results using the large- A'^c limit from the very beginning. At large 
Nc we can apply the saddle point method to the calculation of the integral on the RHS of Eq. H4.6|l . Indeed, the 



integrand contains the factor [Tr (i?^*'')]^'' which leads to the saddle point equation 

<5i^Tr (i?g*' ) - . 

The original integral runs in Eq. (|4.ti|) over the SU (2) matrices R. But the saddle point method leads to a deformation 
of the "integration contour" so that we must look for solutions R and their variation Su in the complexification of 
SU{2), i.e. in SL{2, C). Taking the variation Su in the form of an infinitesimal SL{2, C) rotation 

i?' = (1 + i(5wV")i? (4.14) 

with arbitrary complex infinitesimal parameters (Joj", we can rewrite the saddle point equation 14.13|l in the form 



(4.13) 



Jw'^Tr (T"i?5" ) = . 



(4.15) 



Thus 

Tr (r^Rg^') ^ . 

We see that 

where c is some complex number. Taking the determinant of this equation, we find 

(det R) (det 5) = . 

Since R belongs to SL{2, C), we have 

det i? = f . 

Hence 



det g = c 



Now we find two saddle points from Eq. (|4.17|) : 



(detg)^/'' 

The calculation of the Jacobian of fluctuations around these saddle points leads to the following general result 
I dRf{R) [Tr {Rg^^)]""^ ^J^, [Tr ■ 

Inserting expression ifi:^ for R± and taking f{R) = ZJj (i?"^), we derive from 
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(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 

(4.21) 
(4.22) 



dRDl^^{R-^)[Tr{Rg'^)] 



2 (det; 



.1/2 



Din 



(detg) 



1/2 



(4.23) 



Combining this with the asymptotic expression H4.8|l for ctn^ and using the property of Wigncr functions ljA3|l . we 
find the saddle point result for the integral (|4.t)|) : 



CTN^I dRDl^^iR-') [Tr(i?g*'-)]'^^ ""^^ ^STTT 



1/4 



(2detg)~^/^i^|;3,3 



(detg) 



1/2 



(4.24) 



We see that we have reproduced the above results H4.1U|I . (|4.12|l for w{g) and arTsJsig)- 



V. GENERATING FUNCTIONAL FOR THE BARYON DISTRIBUTION AMPLITUDE 

A. Baryon distribution amplitude 

The definition of the baryon distribution amplitude uses an auxiliary light-cone vector n. This vector n can be used 
to impose light-cone gauge 



(n • A) = . 



(5.1) 



We are interested in the leading-twist baryon distribution amplitude which is determined in terms of the "good" 
components of the quark field q, 

{n-j)q,f{Xn). (5.2) 
In order to separate the "good" components we introduce Dirac spinors Ug associated with the vector n, 



Us (g) Us = i(n • 7)(1 - 2575) f s = ±i ) , 



(5.3) 



9 



and define the fields 

/oo 
—u_,qefiXn)exp[iXx{nP)]. (5.4) 
-oo 

Now we can define the baryon distribution amphtude as the transition matrix element between the vacuum and the 
baryon B with momentum P: 

*(/lsi)(/2S2)...(/N,Sivj(2;i,2:2, • ■ ■,XnJ 

= Y~fc^C2...CN, (0|Xci/lSl {xi)Xc2f2S2 {X2) ■ ■ ■ XcnJn^S^, {xnJ\B{P)) . (5.5) 

This definition of the distribution amplitude ^ is independent of the vector n. Usually it is convenient to normalize 
n by the condition (nP) — 1 but we do not impose this constraint. Indeed, at large Nc we have P — 0{Nc). We keep 
n fixed at large Nc, therefore {nP) grows as 0{Nc): 

n = 0{N^c) , {nP) = 0{Nc) . (5.6) 

At large Nc the distribution amplitude is concentrated at 

xk ~ l/Nc . (5.7) 

Therefore it is convenient to introduce new variables 

yk = NcXk (5.8) 

which behave as 0{N^). 

B. Generating functional 

Now we define 

/•oo /-oo yOO 

'^{9)^^-^"'^ I dyi dy2... rfyw,5/isi (2/1)5/232 (y2) ■•■ (yjvj 

JO JQ Jo 

The factor Nc ^"^'^ is inserted here in order to compensate the Nc growth of the contribution of the kincmatical factor 
{nP)^"/'^ coming to Eq. (|5.5|) from the product of Nc fields x H5-4|) . 
At large Nc we can write the general decomposition (|1.4|) 

$(.g) = N-^Aig) exp [NcW{g)] [l + O (iV"!)] . (5.10) 

From the definition (|5.9|l of $(5) it is obvious that for any constant A we have 

4>(A5) = A^=a>(g) . (5.11) 

Therefore 

H^(Ag) = M^(.g) + lnA, (5.12) 



A{\g) = A{g) . (5.13) 
Differentiating identity H5.12|l with respect to A, we find 



10 



The distribution amplitude contains the momentum conserving delta function: 

*(/isi)(/2S2)...(/iv,s„j(a;i, 2:2, ■ • ■ , xnJ ^ 6{xi +X2 + ■■■+XN,-1). (5.15) 
Because to the presence of this delta function, the functional <&(<?) has the following property. For functions 

(5.16) 



we have 



-,iyo:P) 



This corresponds to the following singularities of W{g) 



ReW 



^{V0,|3) 



ifyo^l, 
00 if yo 1 • 



-00 if yo 7^ 1 , 
+00 if yo = 1 ■ 



(5.17) 



(5.18) 



VI. EVOLUTION EQUATION 



A. Evolution equation for the baryon distribution amplitude 



The evolution equation for the baryon distribution amplitude at arbitrary Nf. can be easily read from the literature 
starting from the original papers [HIB]' However, one should be careful generalizing the standard equations written 
for the 3-quark distribution amplitude to the Ai"c-quark case. In large-A''^ QCD the baryon distribution amplitude 
^'(a;i, . . . , xn^) depends on on Nc variables Xk [one of them can be eliminated using the conservation of the momentum 
(15.15(1 ], In QCD with Nc colors, the dependence of the leading-twist baryon distribution amplitude on the 
normalization point /i is described by the evolution equation 



jVc + lQ,(/i) 

2Nr TT 



J2 K^j^'^ix,,.. 

l<i<j<Na 



(6.1) 



We work in the leading order in the strong coupling as (fJ-) ■ We have omitted the flavor and spin indices concentrating 
on pair-interaction structure of the evolution kernel. The indices i,j of Kij imply that this operator acts on the 
variables Xi,Xj. For example, the compact equation 



should be expanded as follows 

0(/l^l)(/2^2)(2,^^2,2) 



/l^l/2^2 



4> = K^il^ 



The kernel K is diagonal in both flavor fk and helicity Sk- 

rv'(/lSl)(/2S2)/ / ' \ Sjfl S:f2 ijSi CS2 f>SiS2 ( I l\ 

(/isi)(/2si)V^i'^2,a;i,2:2j — 0^, 0^, 0^,^ 0^^ /V yXi,X2,X-^,X2) ■ 

More information on the evolution kernel and its properties can be found in Appendix IbI 



(6.2) 



(6.3) 



(6.4) 



B. Evolution equation for the generating functional 



Now we want to consider the large- iVc limit. As was explained above, a consistent analysis of the large- A'^c limit is 
possible in terms of the generating functional $((7) H5.9|l . Therefore the first step is to rewrite the evolution equation 
(|6.1() in terms of the generating functional ^{g). Let us contract Eq. H6.1|l with the product of functions gf^s^ [yk) 



d_ 



n / dykgf^s^{yk) 

2Nr TT 



(/iSi)---(/n<:SnJ VAr'"'' Nr 



Nc r 



E n 

l<i<j<Na fc=l 



dykghsAvk) 



^ H/iSi)---(/NcSNc) \ TV 



Nc 



(6.5) 
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Note that the operators Kij are invariant under dilations. Therefore the transition from Xk to yu = N^xu does not 
change the form of these operators. AU the N^^N^ — l)/2 terms of the smn over i < j on the RHS are identical and 
can be written in terms of the generating functional $(g) 



E E 

l<i<j<N^ fmSm 

Nc{N, - 1) 



n / dykdhsAv 



E 



[^12^ j(/isi)...(/«.^ivj I 



- dyi dy2 ^ 5/131(2/1)5/232(2/2) ^4:12^- J-^J-, TrT^^lff) 

2 Jo Jo f 7~^, L ^5/isi(yij 05/2^2(2/2) 

/1/2S1S2 



Now the evolution equation Ht).5|l takes the form 

d ^ . ^ 7V, + la,(^) 
/^^$.(5) = - 

where we use the short notation 



2Nc 2tt 



{g®g)-K 



{g®g)-K 

pOO pOO 



%isi(2/i) %2 82(2/5 



= / dyi dy2 5/isi(2/i).9/282(2/2 

Jo Jo ft 

/1/2S1S2 

/"OO /"OO /"OO /"OO 

= / ^2/1/ ^^2/2 / dy[ / E 

Jo Jo Jo Jo f f c » 

/1/2S1S2 

S 6 

X 5/isi (2/1)5/232 (2/2)^^''"' (2/1, 2/2; 2/i, 2/2)7 Tm i rTT*M(5) 

SghsAVi) 05/282(2/2) 



r<&M(.9) 



Here we used Eq. (|6.4|l . 



(6.6) 



(6.7) 



(6.8) 



C. Evolution equation in the leading order of the large- A^c limit 



Now we can study the large- A^c limit of the evolution equation (|6.7|l . We insert the large- iVc ansatz (|5.1(J|I into Eq. 
(|6.7|l and obtain in the leading order of the 1/A^c expansion: 



d_ 



ai^^) \ (9 9) ■ K ■ 



5W,,{g) ^ SW^ 



59 



where 



This limit exists since 



a(^) = lim 



as(p)Nc 
An 



(6.9) 



(6.10) 



(6.11) 



Eq. H6.9(l is obviously compatible with the constraint H5.12|l on 14^^(3). The evolution equation (|6.9(l is nonlinear in 
Wp,{g). Being a first order differential equation in /u, it allows us (in principle) to find at any normalization point 
IJL starting from the initial value for W^^ at some point /zq. 

Instead of jjL it is convenient to introduce the new variable t such that 



dt — 2a (/i) 



(6.12) 



Then 



1^(5,0 



{{g(®g)-K 



m{g,t) ^ 5W{g,t) 



S9 



(6.13) 
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Note that the evolution equation (|6.13() does not depend on the quantum numbers of the baryon. This agrees with 

the statement that functional W^{g) is universal for all low- lying baryons. The difference between baryons appears 
in the pre-exponential factor A{g) H5.10|l . 

D. Hamilton- Jacobi structure 

Eq. (|6.13|l has the form of a functional Hamilton-Jacobi equation. Changing the notation 

g^q, W{g,t)^S{q,t) , (6.14) 

we easily see that we deal with the equation 

where the Hamiltonian (written in terms of discrete variables p„, 

H{p,q) = ^ X! ^ijrnnqiqjPrnPn (6.16) 



2 

ijmn 



is quadratic both in coordinates qn and in momenta 



Pn = 1^ • (6.17) 

oqn 



E. Evolution equation in the next-to-leading order 



Inserting ansatz H5.1()|l into the evolution equation H6.7|l . we can also derive an equation for the pre-exponential 
factor A{g,t) 



d_ 

dt 



\nA[g,t) 



[g®g)-K- 



5W{g,t) 5\nA(g,t) 



+ m {{g®g)-K 



5g Sg 
SW{g,t) ^ SW{g,t) 



where 



b{t) 



2 1(5® 5)-^ 



- Hm < Nr 

2 JVe^oo 



^<^^]Wig,t) 



Sg 

Sg Sg 
Nc + 1 asip) 



47r a(/i) 



(6.18) 



(6.19) 



Taking the difference of two equations H6.18|l for baryons Bi, B2 described by the pre-exponential factors ABi{g,t) 
and AB2{g,t), we obtain the following equation for the evolution of ABi{g,t)/AB2{g,t) 



d Ab, (ff, t) 
dt ABA9,t) 



{g (g>g) ■ K ■ 



SW{g,t) ^ s ABAg.t) 



Sg 



Sg Ab2 {g,t) 



(6.20) 



This is a linear homogeneous equation for Ab^ (5, t) /Ab2 {g,t)- The structure of this equation agrees with the statement 
about the factorization of AB{g,t) into "elementary functionals" ^i{g,t) [cf. Eq. (|3.8fl ]: 



AB{g,t) ^ Aoig,t)Y[[C,ig,t)Y 



(6.21) 



Functionals ^i{g, t) correspond to the elementary excitations Ui in the mass formula H3.3|) . The factorization (|6.21|) is 
compatible with the evolution equation H6.20|l if functionals £,i{g,t) obey Eq. H6.20|l : 



d 

QlU{g.t) = -\{g®g)-K' 



sw{g,t) ^ 5Ug,t) 



Sg Sg 

Strictly speaking, expression H6.21|l must be modified by the zero mode factors which are discussed in Sec. IX CI 



(6.22) 
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VII. ASYMPTOTIC LIMIT 

A. Double limit Nc — > oo, t ^ oo 

At finite Nc the asymptotic large- /j, behavior of the baryon distribution amphtude is well known [see Eq. (|8.3()(l ]. 
In terms of the variable t (|().12|l the limit of large scales /i corresponds to t ^ oo. Now we want to study the large-f 
behavior of the functional W{g, t). Since the functional W{g, t) is defined via the large- iVc limit, we face the problem 
of the double limit N^, ^ oo, t ^ oo. As we shall see, the order in which these two limits are taken is important. 
One of the results of our analysis of the double limit Nc ^ oo, t oo will be a general expression for the anomalous 
dimensions of the leading-twist baryon operators in the large Nc limit. Later (in Sec. IIX|I this expression will be used 
for the calculation of these anomalous dimensions. 



B. Linearized theory of the asymptotic limit 

Now we turn to the investigation of the asymptotic limit of W{g,t) a.t t ^ oo. In this section the analysis will be 
based on the evolution equation H6.13|l . A complete systematic analysis of the nonlinear evolution equation H6.13|l in 
the asymptotic regime is beyond the scope of this paper. Our aim is more pragmatic: we want to find the physical 
solution. On the way to this physical solution we shall use various assumptions about the structure of the solution. 
The justification of these intermediate assumptions comes from the final form of the solution and from the independent 
analysis of the asymptotic limit in Sec. IVIII CI 

We expect the following structure of the large-i behavior of the functional W{g,t) 



W{g,t) W^{g,t) = Woig) - , 

where the function a{t) is g independent and the functional Wo{g) does not depend on t. 

Inserting the ansatz (|7.1|) into the evolution equation (|6.13f) , we arrive at the system of equations 



{g(g)g)-K 



SWoig) ^ SWoig) 



(7.1) 



(7.2) 



a{t) = Et, (7.3) 

where E is some constant. 

According to Eqs. (|5.12|) and H7.1|l the functional Wo{g) must obey the condition 

WoiXg)^Wo{g)+lnX. (7.4) 

Note that in terms of the Hamilton- Jacobi interpretation H6.15|l of the evolution equation (|6.13|) the asymptotic 
equations (|7.2f) and (|7.3|) correspond to fixing the energy E in the Hamilton-Jacobi equation: 

This "energy" interpretation of the asymptotic regime is quite similar to the standard analysis of the asymptotic 
distribution amplitudes of hadrons in terms of the lowest eigenstates of the effective Hamiltonian describing the 
evolution. However, our current large- TVc evolution equation is nonlinear and this leads to certain complications. 

The true physical functional W{g,t) is different from the asymptotic solution VFas(5,^) H7.1|l : 

W{g,t)^WU9,t) + AW{g,t). (7.6) 
Both W{g,t) and Was(g,t) satisfy condition H5.12|l . Therefore we must have 

AWi\g,t)^AW{g,t). (7.7) 
At large t the AW{g, t) becomes small and we can linearize the evolution equation H6.13|l in AW{g, t): 

'6WUg,t) ^ SAW{g,ty 



^AW{g,t)^-{{g<E>g)-K. 



Sg Sg 



(7.1 
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According to Eq. (|7.1|l we have 



Therefore 



d_ 
dt 
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{g®g)-K- 



Sg 



SWoig) ^ dAW{g,t) 



Sg 



In order to solve this equation, let us first consider the spectral problem 

'6Wo{g) ^ 6U9 



{g(g)g)-K- 



5g 



in the class of functionals S,k{g) obeying condition H7.7|l 

U^g)^Ug)- 

The solution of Eq. H7.10|l can be represented in the form 

k 

Now we insert Eqs. (EU, lO, and (TTT^ into Eq. (|T^ : 

W{g, t) '=°° Wo{g) -Et + Y, afce-"''*a(g) , 



(7.9) 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



(7.14) 



This equation is derived in the double limit Nc — > 00, t ^ cxd. Note that the limit Nc — *■ c» is taken first [remember 
that the functional W{g, t) is defined by Eq. H5.1()|l via the large- A^c limit]. 

C. Asymptotic behavior in terms of anomalous dimensions 

Formally Eq. H7.14|l solves the problem of the asymptotic limit t — > 00. Nevertheless one can wonder how this 
solution is connected with the standard representation of the t dependence in terms of the decomposition in operators 
with given anomalous dimensions. For the generating functional ^{g,t) (|5.9() this decomposition has the form 



Hg,t) 



J2^»ig)exp{~rat) . 



(7.15) 



Here Fq, are anomalous dimensions of operators diagonalizing the evolution. The expansion (|7.15l) is written for finite 
Nc- If we take the large- A'c limit in Eq. (|7.15|) . then we must be careful about the order of the two limits A^c ^ 00 
and t 00. The approach based on Eq. H7.15|l assumes that the limit of large t is taken before the limit Nc ^ 00. 
On the contrary, Eq. (|7.14|) deals with the opposite order of the two limits (first Nc ^ 00 and then t 00). 
At large Nc we expect the following structure of the anomalous dimensions Fq, 

F„ = ciA^, + co{a) + c^i{a)N-^ + 0{N;^) . (7.16) 

Note that the leading order coefRcicnt ci is independent of a. This independence will be explained later. 
If we insert decomposition H7.16|) into Eq. (|7.15|) . then the 1/Nc expansion for F^ will be exponentiated: 

$(5,t) *^°°^^^°°V<f„(g)exp{-t[ciA^,-Fco(a)+c_i(a)A^-i + ...]} . (7.17) 



This exponentiation of the 1/A^c expansion for Fq, means that we have a rather nontrivial sensitivity to the order of 
the limits i — > 00 and Nc 00. 
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If one is interested in the other order of hmits (first Nc oo and then t oo), then one should start from Eq. 
(|7.14|) . Inserting the decomposition 1)7.1411 into Eq. (|5.1U|I . we find 



Woig)-Et + Y, 



ate 



(7.18) 



Now we can compare the two asymptotic expressions (|7.17|l . (|7.18|l derived for functional $((?, i) in the double limit 
Nc — !■ cx), i — !■ cx) assuming different orders of these two limits. Separating the factors containing the exponentiated 
product Net, we conclude from the comparison of Eqs. (|7.17|) . H7.18|l that 

(7.19) 
(7.20) 



We see that 



ci=E. 



This result explains why the coefficient ci appearing in large- A'c expansion (|7.1t)|l for Ta is independent of a. 

One can also match the two asymptotic representations (|7.17|l and (|7.18|l in the next order of the 1 /N^ expansion. 
To this aim we expand the RHS of Eq. H7.18|l in powers of afce~^'=*^fe(g): 



Hg,t) ^'-^'--^ A{g,t)e^v{NcWo{g)} 



X ^ exp 



Comparing this decomposition with Eq. H7.17|l . we see that 

co(a) = A£; + ^nfcf7fe (nfe = 0, 1, 2, . . .) 



(7.21) 



(7.22) 



where the 0{Nc) contribution AE comes from exponential part of the large-t asymptotic behavior of A{g, t) ^ e 
Inserting this expression for co(q;) into Eq. H7.16|l and using the expression H7.20|l for ci, we find 



r{„,} = NaE + (^AE + rik^ 



k +c,iinun2,...)N-^ +0{N-^). 



Here we label the anomalous dimensions Ta with sets of integer numbers {rife}: 

a = {nk} (nfe = 0,1,2,...). 



(7.23) 



(7.24) 



At finite Nc the spectrum of the anomalous dimensions comes from the diagonalization of the A^c-particle "effective 
Hamiltonian" and the anomalous dimensions Fq, are labeled by Nc — 1 "quantum numbers" a 1,0:2, . . .aN^-i (one 
degree of freedom is eliminated by the momentum conservation). In the limit Nc — > 00 we arrive at the "oscillator 
spectrum" ()7.23|l and the anomalous dimensions are parametrized by a set of integer excitation numbers Uk- 

The dependence of the anomalous dimensions on some quantum numbers (e.g. helicity) appears only in the order 
0{Nc^)- Denoting these quantum numbers by p, we can write the correct generalization of Eq. H7.23|l 



r{nfc,p} - NcE 



AE + Y, rik^k + c_i(ni, 712, . . . ; p)N-^ + 0{N-^) . 



(7.25) 



This structure of the spectrum corresponding to an effective harmonic oscillator in the order (with anharmonic 
corrections appearing in the next orders of the 1/Nc expansion) is typical for semiclassical systems. The semiclassical 
nature of the \/Nc expansion is known since long ago. Depending on the methods used for the construction of the 
1 /Nc expansion, the semiclassical features can arise via the mean field approximation, saddle point approximation in 
path integrals, WKB expansion etc. 

Our derivation of the expansion (|7.25() was based on the asymptotic linearized analysis of the Hamilton- Jacobi 
equation H7.10|l in terms of the large-^ decomposition H7.14|l . In this approach the spectrum of Vlk is given by the 
eigenvalues of the spectral problem (|7.11|) for the functionals £,k{g)- 

In the rest of the paper we compute parameters E, AE and f2fc appearing in the expansion H7.25I) . The results are 
given by Eq. H8.9|l for E and Eq. H9.4|l for AE. The parameters are computed in Eqs. H9.40I) . H9.41|l in the case 
of one quark flavor Nf = 1. The role of parameter p is played by the helicity J3 in this case. The 0{N~^) term 
c_i(0, 0, . . . ; Jz)Nc^ is computed for the lowest anomalous dimension corresponding to the asymptotic wave function 
of baryons with given helicity J3 [see Eq. (|9.2|) ] . 



16 



VIII. CALCULATION OF THE ASYMPTOTIC FUNCTIONAL ^^0(5) 

A. Methods 

We want to compute the asymptotic functional Wo{g) which determines the asymptotic behavior 1)7. l|l of W{g,t) 
at large t. Two methods can be used for the calculation: 

1) We can solve Eq. (fT^ . 

2) We can use the traditional approach to the determination of the large-t asymptotic baryon distribution amplitude 
at finite Nc- Using this finite- A'^c asymptotic distribution amplitude we can construct the corresponding generating 
functional by analogy with Eq. I|5.9|l . The large- iVc asymptotic behavior of this functional will be dominated by the 
exponential term containing information on Wo{g). 

The first method is described in Sec. IVTTT Rl In Sec. IVTTTd we show how the second method works. In Sec. IVIII Dl 
we compare these two methods. The computed functional Wo (5) depends on the number of quark flavors Nf. The 
nature of this dependence is studied in Sec. IVIII El In Sec. IVIII Fl we comment on the analytical properties of the 
fmictional W{g,t) and illustrate the general statements using our result for the asymptotic functional WQ{g). 

B. Solution of the equation for Wo{g) 

According to Eq. IjT.ll) the main quantity characterizing the asymptotic limit is the functional Wo{g). This 
functional can be found by solving Eq. (|7.2() . 

Our approach to Eq. (|7.2(l will be rather heuristic. We simply use an ansatz which allows us to solve Eq. 1)7. 2|l for 
some special value of E. The relevance of this value of E for the description of the asymptotic behavior will be seen 
from the analysis of Sec. IVIII CI The form of the solution Wo{g) of Eq. H7.2|l depends on the number of quark flavors 
Nf . We start from the simplest case Nf = 1 . The result found for Nf = 1 is easily generalized to the more interesting 
case Nf = 2. 

1. Case Nf ^1 

In order to solve Eq. I|7.2(l we use the property of the evolution kernel K^'^^^ 



which is valid for any functional X{g). This identity follows from Eq. I|B21|1 derived in Appendix^ Keeping in mind 
this identity, we make the following ansatz 

'^=ye--i^^yU.(g), (8.2) 
09s{y) 

where X{g) and Us{g) are some functionals of gg. Inserting this ansatz into Eq. (|7.2() and using Eq. H8.1|l . we find 

Y,[Gs{g)Us{gf ^AE, (8.3) 



where 



poo 

Gs{g)^ / <?.(y)ye-^(^)^. (8.4) 
"'0 



Combining this expression for Gs{g) with Eq. (|8.2() . we find 







' dygs{y)^-pl^=Gs{g)Us{g). (8.5) 



According to Eq. I|5.14|l we have 
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Combining the last two equations, we conclude that 

s 

We can satisfy both Eq. 1)8. 3|l and Eq. 1)8. 7|l by taking 



which corresponds to 



Now it follows from Eqs. (IH3, 1H3I), and 



^^0(9) ^ i-Xig)y 

Sgsiy) 2^ 



E 



dy'gs{y')y'e 



We can look for the solution of this equation in the form 

1 r f°° 

Woig) = F [Xig)] + 3 ^ l^^i, "^y^^^^^^^"''^'^ J ' 

where F{X) is an arbitrary function and the functional X{g) is imphcitly determined by the equation 
o I 1 r poo 1 I 



(8.7) 



(8.8) 



(8.9) 



(8.10) 



(8.11) 



.12) 



Thus we have found an infinite set of solutions Wo (.9) of Eq. (|7.2(l corresponding to the same constant E ()8.9(l . These 
solutions are parametrized by arbitrary functions F{X). The physical solution is fixed by the condition H5.18|l which 
leads to the following choice of F{X): 

F{X) = X. (8.13) 
Indeed, taking this function F{X) and computing the LHS of Eq. H8.12|l for functions g^yo^l^^ ()5.1t)f) . we obtain 



f) ( 1 r r°° 



= 1 - yo ■ 



.14) 



Therefore Eq. H8.12|l has no solutions at yo 7^ 1 in agreement with the condition 15.18|l . 
We can combine Eqs. (|8.11|) . H8.12|l . and (|8.13|l into the final representation for Wo{g) 

{1 r f°° 

X + 7;y2^^ / d,ygs{y)ye~^'" 
^ , Uo 

The word extremum must be understood in the sense of the equation 

1 r f°° 



.15) 



d 

dX 



-Xy 



= 0. 



(8.16) 



In the general case both gs and X can be complex so that we cannot speak about a maximum or a minimum. Equation 
(|8.16l) determines the functional X{g) and this functional should be used on the RHS of Eq. H8.15() for the calculation 
of 1^0(5)- 
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2. Case Nf = 2 



The generalization of the above results to the case Nf = 2 is straightforward. By analogy with Eq. (|8.2fl we make 
the ansatz 



Introducing the notation 



we find the Nf = 2 analogs of Eqs. (IH31) and 



One can satisfy both equations by taking 



(8.17) 

(8.18) 

(8.19) 
(8.20) 

(8.21) 



which leads to the same value 



E 



as in the N f ~ 1 case (|8.9f) . 

According to Eqs. (|8.18() and H8.21|l we have 

Ufs{g)^\[G{g)]:} = \ 

Inserting this result into Eq. H8.17|l . we find 

— = -ye ■ 

ogfs[v) 2 



dy9{y)ye 



-X{g)y 



dy'g{y')ye-''^^~>y' 



The general solution of this equation compatible with the momentum conservation constraint H5.18|l is 

dygfs{y)ye~'^'y 



Wo(q) = extremum < X ^ — Indet 
X I 2 fs 



.22) 



.23) 



(8.24) 



(8.25) 



Similarly to the Nf = 1 case, this compact representation should be understood in the sense of the (generally speaking 
complex) "extremum equation" 



d 

dX 











In det 


[/ 




fs 


Jo 



dygfs{y)ye 



0, 



which determines the functional X(g). Eq. H8.26|l can be rewritten in the form 

n -1 ■ 



Tr 



dyg{y)y e 



dy'g{y')y'e-''y' 



2. 



Inserting Eq. (|8.25(l into Eq. 1)7. 1|1 . we find the final expression for the asymptotic functional: 

1 



Was (5, i) = extremum < X + - In dot 
X I 2 fs 



dygfs{y)ye 



Et. 



(8.26) 



.27) 



(8.28) 
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C. Direct construction of the asymptotic functional 

Let us consider the case of Nf — 2 flavors. The asymptotic distribution amphtude of baryons with heUcity 

I J3I < T 



(8.29) 



c, 



isi)(/2S2)...(/«.s„j2:ia:2 . . . xnA^i + + . . . + xn^ - 1) . 



(8.30) 



The hehcity-flavor part ^'l^^si){f2S2) (/» ) '^oi^'^i'^'^s with the wave function of the toy quark model H4.2|) up to a 



normalization factor: 



C, 



(/iSi)(/2S2)---(/NcSNc) 



(/lSl)(/2S2)...(/jVeSNc) 



.31) 



3A^ /2 . . . 

The factor of Nc " is inserted in order to simplify the subsequent expressions (the question about the large- A'c 
behavior of the coefficient a!--^^ requires a separate analysis). 
Combining Eqs. (|121), (jHilB, and (pOT|l . we find 



(B) 

(JlSl)(j2S2 



k=l 



X1X2 ■ ■ ■ XN^ixi + X2 + ■ ■ ■ + Xn, - 1) • 



.32) 



It is easy to see that this wave function diagonalizes the RHS of the evolution equation (|6.1|l . Indeed, using identity 
(|B20|) . we find 



^*J*(/lSl)(/2S2)...(/«,5ivJ - 2 '''"^ (/l«l)(/2«2)...(/iV,S„J 



(8.33) 



The helicity-flavor wave function C'^^/'/sj)(/2S2) (/„ sjv ) ^^"^ baryon with helicity J3 has + J3 indices Sfc = +1/2 
and — J3 indices Sk — —1/2. Therefore 



L)(/2S2)...(/iV,S„J - 2 (/l"! 

l<i<J<Wc l<i<j<A'e 



)(/2S2)...(/NcSNc) 



(f + J3)(f + J3-I) , (f -J3)(f -J3-I 



(/iSi)(/2S2)---(/NcSNc) 



fN, 



- 1 + ( J3) 



(/iSi)(/2S2)...(/NcS1Vc) 



Therefore the evolution equation H().l|) for the wave function (|8.3U|) becomes 



1 



iJ3 



'^''{xi,...,xnJ. 



Taking the limit of large Nc and assuming that J3 = 0{N^) in this limit, wc find 



,xnJ = -Y^as{iJ-)'i/^{xi,. . . ,xnJ ■ 



This corresponds to the fj, dependence 

•^^'{xl,...,XNJ^ (a;i, . . . , ijvj exp {Nc [a{^io) ~ cr(M)]} 



where function cr(/i) obeys the equation 



da{ii) ^ Nc_ 
167r 



.34) 



.35) 



(8.36) 



..37) 



(8.38) 
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Changing from /i to the variable t (|6.12() we find 



a = -t. 



This agrees with the above results (|7.3|l and (|8.22|) . 

Thus the evolution of the asymptotic distribution amplitude (|8.32|l is described by the /i dependent factor 



so that 



(/lSl)(/2S2).--(/]V<:SiVe) 



(xi, X2, . . . , xnJ = a^^^CTN.N^^-/^ exp [-N,a{n)] 



(8.39) 



(8.40) 



XiX2 ■ ■ ■ Xn^5{xi + X2 + ■ ■ ■ + xn^ — 1) 

dZ 



2tt 



..41) 



fc=i 



Next we compute the generating functional (|5.9() for this wave function 



$.(.9) ';^/-^- exp[-iV.a(^)] JdRj dZDl^.^{R-^) 



X exp <^ Nc < -iZ + In 



Rfs / dygfs{y)ye'^'^ 



..42) 



The integrals over R and Z can be taken using the saddle point method. Let us first perform the integration over R 
at fixed Z. This can be done by repeating the calculation of Sec. IIV CI for the quark model. We must simply replace 



5 ^ / dyg{y)ye'^y 



in Eqs. (|4.6() and (I4.21() . As a result, we find the saddle point for the R integral in Eq. 18.42|l : 



dy9{y)ye 



iZy 



det 



dyg{y)ye 



iZy 



-1/2 



(8.43) 



(8.44) 



Inserting this result into the exponent of (|8.42|l . we find with the exponential accuracy (neglecting factors independent 
of g and /i) 



$-(.g) ^ exp [~NMp)] 



X I dZ exp I -iZ +2^^ 



dyg(y)ye'^y 



The saddle point equation for the Z integral is 



d_ 

dZ 



—iZ + - In dot 



dyg{y)ye'^y 



= 0. 



This equation coincides with Eq. (|8.26|) if we set 

X = -iZ. 

Now we see that the saddle point integration in Eq. 18.4511 reproduces our old results H8.25|l and 



.45) 



(8.46) 



..47) 
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D. Comparison of the two methods 

Thus the two methods described in Sections lVIII Bl and lVIII Cl have led us to the same resuh (|8.25|) for the functional 
Wo (5). Note that in the second method we assumed that the baryon obeys the condition | J3I < T which was needed 
in order to make use of the nonrelativistic spin- flavor wave functions H8.31|) . In principle the condition | J3I < T is 
automatically satisfied for the lowest 0{N~^) excited baryons H3.2() with J — T. But the discussion of the universality 
of the functional W{g) in Sec. 11111 shows that the same functional W{g) must describe not only the lowest J = T 
baryons H3.2|l but also higher 0(iV°) excitations and baryon-meson scattering states. Note that the first method 
based on the asymptotic analysis of the evolution equation for W{g) is completely compatible with the universality of 
W{g). Actually the condition | J3I < T used in the second method can be omitted. Indeed, the origin of this condition 
is the totally symmetric form of the Xk dependent part of the asymptotic wave function H8.30|l . In the case of the 
baryons with | Jaj > T we cannot keep this form of the Xk dependence. But one can construct the corresponding 
asymptotic wave function and show that it is "almost symmetric" in the sense that in the large- A^c limit this wave 
function leads to the generating functional which has the same W{g) exponent, and all effects of the asymmetry 
of the wave function are localized in the pre-exponential factor A(g). 

One of the advantages of the second method considered in Sec. IVIII CI is that it allows us to trace the connection 
with the traditional analysis of the asymptotic baryon distribution amplitude. Another good feature of this method 
is that it clarifies the saddle point origin (|8.44|l . H8.46|l of Eqs. H8.24|l . H8.26|l which appear rather formally in the 
first method. Concerning this saddle point interpretation, we must make an important comment. In principle the 
large- A^c limit justifies the applicability of the saddle point method to the calculation of the integral H8.42|l . However, 
the large-A^c limit does not guarantee that 

1) the saddle point equation H8.46|l has one and only one solution, 

2) the R and Z "integration contours" can be properly deformed so that the solution of the saddle point equation 
gives the dominant contribution. 

Generally speaking, the validity of the conditions needed for the applicability of the saddle point approximation 
depends on the function g for which we compute the asymptotic functional Wo{g). It is easy to construct examples 
of functions g for which the saddle point method works and counter-examples illustrating the violation of the saddle 
point method. 

To summarize, our representation for the asymptotic functional (|8.28() should be considered as a formal expression 
that carries information about the saddle point equation but does not fix the relevant solution of this equation. 



E. Nf dependence 



The results of our analysis of the cases Nf — 1 and Nf = 2 are represented by expressions (|8.15(l and (|8.25(l for the 
functional Wo{g). We stress that the Nf — 1 and Nf — 2 results are different. Let us use the notation W^^^\g^'^f^) 
in order to mark the Nf dependence. Starting from some gi^\ let us define 

One could naively expect that this choice of gj:J selects only the contribution of the quarks with / = 1 so that the 

values VFo^"* (g'^') and W^^^g^^^) must coincide. But this does not happen. Moreover, the functional Vt^o^''((7'^') 
has a singularity for the functions (|8.48|) because of the vanishing determinant 

det I dygf}[y)ye-''y 



(8.49) 



ReWo^^^(5(2)) = -cx). (8.50) 

What stands behind this difference of the N f = 1 and Nf = 2 cases? Obviously the Nf dependence of the nonpertur- 
bative dynamics of QCD has nothing to do with this problem because the evolution equation does not know anything 
about this nonperturbative dynamics. This difference has another origin. 

The large- A'c baryons which we studied in the Nf = 1 case consist of quarks with the same flavor u. If we "embed" 
this {uu . . .u) baryon into the Nf = 2 theory, then it will be classified as a state with isospin T3 = Nc/2. Thus our 
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large- A^c work in the Nf = 1 sector should be interpreted in the Nf = 2 terms as the limit 

^ oo , Tg = ^ ^ OO . (8.51) 

However, the functional Wg^^^ls^^'') introduced for another limit: 

Nc^oo, T = 0(A^°) . (8.52) 

Although we did not emphasize the condition T — 0{N^), it was really essential in our Nf = 2 analysis. For example, 
in the saddle point calculation of the R integral (|8.42(l we did not include the function D'j^rp^{R'^) into the saddle 
point equation, which is allowed only if T is kept fixed in the limit — > oo. 

Now we understand the meaning of the singularity H8.50|l . After the exponentiation this singularity leads to zero: 

$(2)(g(2))^^(2)(g(2))expJ^^^^(2)(^(2))j (8 53) 

The origin of this zero is obvious: function g^^^ (|8.48l) selects the states with T3 — Nc/2 whereas the baryon has a 
finite isospin which is kept fixed at large Nc- This discrepancy leads to the vanishing functional <I>*^2)(|(j,(2)) 

Thus our analysis of the two cases Nf = 1 and Nf = 2 used different assumptions about the behavior of the baryon 
isospin at large Nc- 

Now we can turn to the case of Nf = 3 flavors. Our experience with the isospin shows that one has to distinguish 
between two cases: 

1) fixed strangeness S — 0{Nc), 

2) growing strangeness S — 0{Nc). 

The case of fixed strangeness can be trivially reduced to the Nf — 2 case. This follows from the universality of 
the generating functional W{g). Indeed, in the large-A^c world with Nf — 3, the nonstrange baryons and baryons 
with 0{Nc) strangeness are described by the same functional W^^^g^^"^). For the nonstrange baryons the functional 
(j)(3)(g(3)-) ||5 9|| obviously does not depend on the strange components of g'^f^ with / = 3. Now we apply Eq. (|5.10() 

to these gg^' -independent functionals $(^^(g('^^) and conclude that the functional W^'^^ {g'^'^^) is also (/g^'' -independent. 
Due to the universality this property of W'-^'' {g'-^'') can be extended to the case of baryons with 0(iV") strangeness. 

Thus in the Nf — 3 world, the baryons with 0(iV°) strangeness are described by the functional W''^\g''^^) which 

(3) (3) 
is independent of the components 53^ . The dependence on g^J appears only in the pre-exponential functionals 



F. Analytical properties of W{g) 

According to Eq. (|5.9|) ^{g) is a homogeneous polynomial functional of g. In particular, <&((?) is a holomorphic 
functional of g. Does this mean that the functional W{g) defined by Eq. H5.10() is also holomorphic in 17? The situation 
is rather subtle. The expression for W{g) via $(5) following from Eq. H5.10|l . 

W{g) ^ Imi -r^ln$(5), (8.54) 

contains a logarithm which can lead to a violation of the analyticity. In principle, the appearance of this logarithm is 
a mere notational effect: instead of Eq. (|5.10() we could work with the large- TVc asymptotic representation 

$b(.9, Nc) = KAsig) [Eig)f^ [l + 0{N-')] . (8.55) 

Obviously 

E{g) = _ (8.56) 

Using Eq. H8.55|l , we find 



(8.57) 
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Here Nc is shifted by two in the numerator of the RHS. This is necessary for the cancellation of the factors Asig) 
Indeed, according to Eq. (|3.1|) the quantum numbers of baryons B are different in theories with odd and even iVj, 
Eq. H8.57|l shows that the functional [E{g)]^ is well defined and has no ambiguities in contrast to the functional W{g) 
In order to illustrate these problems let us consider the toy quark model described in Sec. IIVI Our result H4.10|l for 
the toy analog 'w(g) of the functional W{g) contains a logarithmic singularity for degenerate matrices g with detg ~ 
According to Eq. (|4.1UI) in the toy model we have the analog expw{g) of the functional 



e{g) = exp wig) = ^2 detg . (8.58) 

We see that the analyticity of e{g) is broken by the square root singularity that appears for degenerate matrices g. As 
was mentioned in Sec. IIVI the sign uncertainties of this square root are insignificant for [e{g)]'^'' at even Nc. For odd 
Nc this sign ambiguity is compensated by a similar root ambiguity of the pre-exponential D function in Eq. (|4.12|) . 

Now let us turn to the analytical properties of the functional Wo{g) describing the asymptotic behavior of the 
baryon distribution amplitude at high normalization points fi — > oo. To be specific, let us consider the case Nf = 2 
when Woig) is given by Eq. H8.25|l . It is easy to see that this expression for Wo{g) has two ambiguities. 

1) Expression (|8.25|) contains the logarithmic term 



- In det 

2 fs 



dygfs{y)ye 



(8.59) 



which has an additive 7rm ambiguity. This logarithm is accompanied by the coefficient 1/2 so that the uncertainties 
of exp [A^c 1^0(5)] reduce to (±1)^"^ and disappear in $5(5) in the same way like in the case of the toy quark model. 

2) Parameter X in Eq. (|8.25l) is defined implicitly via the extremum equation (|8.26|l . The ambiguities of the 
determination of X from this extremum equation were discussed in Sec. IVIII Dl These ambiguities can also lead to 
a violation of the analyticity of the functional Wo{g)- 

Note that the appearance of singularities in W^ig) and the violation of analyticity are natural from the point of 
view of the WKB method used in our large- A^c analysis. Such well-known "singular" features of the WKB method 
like caustics. Stokes lines, etc. have their analogs in our problem of the baryon wave function at large Nc- 



IX. DIAGONALIZATION OF ANOMALOUS DIMENSIONS 

A. Anomalous dimensions of baryon operators at large Nc 

We already know that at large Nc the anomalous dimensions ^{uk} of baryonic operators are given by the relation 
()7.25(1 . Parameter E appearing in this equation is determined by Eq. H8.9|l so that 

r{n.,j3} = \Nc + (^^E + Y,nk^})j +c_i(ni,n2,...; J3)^-^+0(7V-2). (g^^^ 



Here J3 is the helicity of the corresponding operator. In the case of one quark flavor N f — 1, the helicity J3 is the 
only quantum number which does not contribute in the order 0{N^). 

The lowest anomalous dimension corresponds to the case when all numbers Uk vanish. We know the corresponding 
anomalous dimension exactly from the evolution equation for the exact asymptotic wave function (|8.35|l : 

r,„..„,.l(^_.),^. 



Note that we do not include the factor 

Nc + 1 ajfi) 
Nr 2tt 



(9.3) 



of the evolution equation H().7|l into the anomalous dimensions (19. 1)9. 2|l and treat this factor separately. The exact 
result for the anomalous dimension r{^„j,=o} l|9.2(l agrees with the general structure of the expansion 1)9. 

Comparing the exact result 1)9. 2|l with the 1 /Nc expansion H9.1|) , we can find the parameter 



(9.4) 
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B. Traditional approach to the calculation of anomalous dimensions 

The "excitation energies" ftk appearing in the large- iVc representation H9.1() for the anomalous dimensions of the 
baryon operators can be computed by solving Eq. (|7.11|l . Our derivation of this equation was based on the linearized 
analysis of the asymptotic behavior H7.14|l of the solution W{g,t) to the evolution equation (|6.13|l . Let us show how 
the same result comes from the traditional approach to the calculation of anomalous dimensions. 

In the traditional approach one has to diagonalize the effective Hamiltonian corresponding to the evolution equation 
(EH) 

H^a{xi, ■ ■ ■,XnJ = Ta'^aixi, ■ ■ ■,XnJ , (9.5) 



H 



l<i<j<Na 



(9.6) 



Let us show that the eigenvalues Fq, appearing here coincide with the expression (|9.1|l at large N^. First we rewrite 
Eq. H9.5|l in terms of the generating functional for the function ^a{xi, . . . , xat^). By analogy with Eq. H5.9|) we define 



Jo Jo 

X *a,(/iSi)(/2S2)...(/iV<,s«J 



We 1(2/1)5/2^2(^2) ■ • ■5/iVeS»e(yiVc) 





yi_ y2_ 



Repeating the same steps that were used in the derivation of Eq. H6.7|l . we find 



2 (5 <^ 5) • ^ ■ 



^(^^ ) $a(3) 



Fa$a(g). 



Next we write the large- A^c ansatz H5.10II 

$„(g) = K-A^ig) exp [N.Woig)] [l + O (N-')] 

and expand 

T^^N,E + AT^ + 0{N-^) . 



(9.7) 



(9.8) 



(9.9) 



(9.10) 



We use the tilded notation ^'q, $0(5), Aa{g) in order to distinguish the wave functions and functionals associated 
with the diagonalization problem H9.5|) from the wave functions and functionals $(5), ■A{g) corresponding to the 
physical baryons states. However, note that the functionals Wo{g), ^k{g) appearing in this section are the same as 
before. 

Inserting ansatz l|9.9|l into Eq. 1)9. 8|l . we reproduce Eq. (|7.2|) in the leading order of the l/Nc expansion. In the 
next order we obtain the equation 



1 



{g(g)g)-K' 



[g(dg)-K- 



5Wo{g) ^ SlnA^ig) 



— (8) — 
<5g Sg 



Wo{g) 



AF. 



In addition to this equation the functionals Aa{g) must be homogeneous, 

AaiXg) = ia(.g) , 
for any number A. This condition follows from the property 

l>„(A5) =A^^I>„(5) 

and from Eq. (|7.4|l . 



(9.11) 

(9.12) 
(9.13) 
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Let us use the value a = for the ground state. Taking the difFerence of the equations (|9.11ll with a ^ and with 
a = 0, we find 



^9 Sg Aoig) 



Now let us make the ansatz 



k 

where ^k{g) are solutions of Eq. (|7.1H) . Then 
According to Eq. (|7.11|) we have 



Inserting Eq. H9.17|l into Eq. ()9.1fci|) . we find 



og og 



= Ar, - AEo . 



Now we combine this with Eq. ()9.in|l : 



AT a - AFo = ^ Uk^k 



Taking a = {nfe}, we see that this result agrees with Eq. (|9.1(l and 

AFo = AE . 



(9.14) 



(9.15) 



(9.16) 



(9.17) 



(9.18) 



(9.19) 



(9.20) 



Thus at large Nc the diagonalization of the effective Hamiltonian H9.5() reproduces our old expression for the anomalous 
dimensions of baryon operators (|9.1|l . 



C. Calculation of excitation energies Q,k 

Now we want to find the excitation energies Q,k- This can be done by solving Eq. I|7.11|l . In this paper we solve this 
equation in the simplest case of one flavor Nf = 1. In order to construct this solution we need auxiliary functionals 



MUg)^ / dygs(2/)/+^e-^(s)^ 



(9.21) 



Tks{g) 



Mks{g) 



MM ■ 

We shall see that the solutions (.?) of Eq. (|7.11() can be constructed as polynomials 

U{g)^Pm[T{g)]. 

Obviously we have for any constant A 

TksiXg) = Tks[g) ■ 

Therefore the functionals ^m{g) constructed as polynomials of Tks (|9.23|l will obey condition (|7.12|) . 



(9.22) 

(9.23) 
(9.24) 
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Assuming the structure (|9.23|l of the solution, we find 



so that Eq. H7.11|l can be rewritten in the form 



E 

ks 



dTks 



5W^{g) ^ STksjg) 
Sg Sg 



^n„^P„,[T{g)]. 



(9.25) 



Below we establish an important property of functionals Tks{g)'- 



{g(^g)-K 



Sg Sg 

k 

Tks{g) + E E B'^:Tr.sAg)Tk-nAg) 



Si 71 — 



(9.26) 



where -B^Ji" are some coefficients. The derivation of the equality (|9.26|) and the calculation of the coefficients B^^}^ 
requires some work [see Eq. (jC31ll in Appendix EJ. But once this work is done, the solution of Eq. (|9.25|l becomes 
simple. Inserting Eq. (|9.26|) into Eq. H9.25|l . we find 



4 4" dTks 

ks 



We look for the polynomial solutions H9.23|l 



~Tks + E E ^kli^nsiTk-r 



Si n—0 



^rnPrn (^) ■ 



{k.,s,} 



■Tk„ 



(9.27) 



(9.28) 



Let us assign the "weight" k to the functional Tks- Then the product Tk-^siTk2S2 ■ ■ ■ Tk„s„ will have the weight 

TO fci + fc2 + . . . + fc„ . (9.29) 

Obviously Eq. (|9.27() is closed with respect to polynomials Pm{T) of any given weight to. Therefore we can classify 
the solutions according to their weight (|9.29|) . This is the reason why we identify the notation to for weight 
and the index to of the solution Pm (T) . Let us concentrate on the linear terms Tms in Pm (T) : 



PmiT) — ^ amsTms + (nonlinear tcrms) . 

s 

Separating these linear terms in Eq. (|9.27|) . we find 



(9.30) 



Tms + {PmmPrnsi + BmoPri 



We have taken into account the property 



Tosig) = 1 



(9.31) 



(9.32) 



which is an obvious consequence of Eq. (|9.22() . The coefficients i?^ depend only on the relative orientation of the 
helicities si, S2 so that we can use the notation 



van mn ' 



r) — s,s o — 

mn mn ' 



(9.33) 



Then 



^ ^ {Bmm'Pnisi + B^^Tms) — {Bmm + -^mo) Pnis + (BmmPm.-s + Bj^qT^ 



(9.34) 
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and Eq. (|9.31() takes the form 

s s 

Obviously we have two solutions for a given weight m. One solution is even in s 



ras ^ni, — s ' 



and the other is odd 

According to Eqs. H9.23|l and (|9.8()(l these two possibilities lead to the solutions 

e±(.9) = P^ATig)] = [Trn+{g) ± Trn-ig)] + [tcrms nonlinear in Tfc, (g)] . 
Parameters il^ associated with these solutions can be found from Eq. H9.35|l : 



(9.35) 

(9.36) 
(9.37) 
(9.38) 

(9.39) 



The explicit expressions for can be found in Appendix O [see Eqs. HC36|I - (|C37|l ]. As a result, we 

find 



1 



1 1 



1 



2 m + 2 (m + l)(TO + 2) 



m+l 



(9.40) 



i=2 



0+ = o- 

m[m + 2) 

The corresponding solutions of Eq. H9.25|l will be denoted Cm(5)- 

{g®g)-K- ® =f^m^m(c/)- 

In the above derivation we had m > 1. However, one should keep in mind that 

0+ = Q. 

The solution of Eq. (|9.25f) for this zero mode is 

C+(g)=ri,+ (g)+Ti,_(g). 
Computing the derivative with respect to X in Eq. (|8.16(l . we find 

E^i«(5)=2. 



This means that 



so that Eq. ijO^ for ^j^(g) 



is satisfied in a trivial way. 



{g®g)-K- 



it{9) = 2 



Sg Sg 



= 



(9.41) 

(9.42) 

(9.43) 
(9.44) 

(9.45) 
(9.46) 
(9.47) 
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In Appendix [HI we find another zero mode of Eq. H7.11|l . According to Eq. ljC20|) we have 

'5Wo{g) ^ <5Crot(g)" 



-0, 



where 



6ot(ff) 



Mo+(g) 



(9.48) 



(9.49) 



This zero mode corresponds to the spin rotations around the space component of the Ught-cone vector used in 
the definition of the distribution amplitude (|5.4|) . (|5.5() . The role of this zero mode will be discussed in Sec. IIXDI 
Let us summarize. The spectrum of the anomalous dimensions of baryon operators is described by Eq. (|9.1|) . 



8 ' 



E 

m>2 



E 

m>l 



+ o(iv-i), 



(9.50) 



with the parameter AiJ ()9.4|1 and with the excitation energies il^ H9.40|l . (|9.4HI . 



D. Helicity 

Our result for the anomalous dimensions (|9.50() is degenerate in helicity J3. This degeneracy is lifted in the order 
0{N~^). It can be seen from the exact (in Nc) result H9.2|l for the operators corresponding to = 0. In fact, the 
conservation of helicity allows us to diagonalize the evolution in terms of generating functionals Aa{g) H9.15|l without 
computing the 0(N~^) corrections to the anomalous dimensions. 

Under the axial rotations corresponding to the helicity 

g's = Qse^'l' , (9.51) 
the functionals M^sig) H9.21|l and T^sig) (|9.22|) transform as follows 

Mks{9')^e''^'M,,{g), (9-52) 
Tksig') = Tksig) ■ (9.53) 

The solutions (,m{g) ^-re constructed as functions of Tks{g) [see Eq. (|9.23() ]. Therefore these solutions are invariant 
under rotations H9.51|l 

^tig') ^ ^taia) ■ (9.54) 

Using Eq. (|9.52|l . we find for the zero mode (|9.49() 

erot(ff') = e'*6ot(5)- (9.55) 

In the previous section we showed that the anomalous dimensions of baryon operators diagonalizing the evolution 
can be parametrized by 

«=(n±,J3), (9.56) 

where are integer numbers and J3 is the helicity of the baryon operator. Assuming this meaning of the subscript 
a in the generating functional (|9.9|l 

A^{g) = A^^^ j^^ig) , (9.57) 

we can write for the rotations H9.51|l 

^{«±,J3}(3') = e'^"^^^{„±,j3}(ff)- (9-58) 
Now let us rewrite the decomposition ()9.15|) in the form 

hnt,.h}i9) = ^(5) I n \ \ n [C(.9)]"" i %ot{g)Y' ■ (9.59) 

Here we have separated the contribution of the zero mode ^rot(5) l|9.49|l . Applying transformation H9.51|l to the 
decomposition H9.59|l and using Eqs. H9.54|) . H9.55|l . and (|9.58f) . we find 

^(5')=^(ff)• (9.60) 
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E. Functional ^0(5) 



Now we want to compute the functional ^0(5) appearing in the decomposition 
for this functional. Taking Eq. (|9.11|) for the ground state a — and using relations 



{g®g)-K 



5W^{g) ^ (51nio(g) 
<5g 5g 



Y^^Y^]W.{g) 



1 

'a 



The solution of this equation can be written in terms of functionals Tks IpH^ : 

M9) = \Y.{T^^^9)-[TU9)f]\ ■ 



First let us derive an equation 
we obtain 



(9.61) 



(9.62) 



A straightforward but tedious calculation [using relations (|9.26ll and (jCSp ] shows that this functional obeys equation 
(|9.61l) . Note that Eq. (|9.61|) does not fix the solution Ao(g) completely. For example, due to Eq. H9.48|l we have the 
freedom of transformations 



io(g) ^ io(g) Krot(5)]' 



(9.63) 



The ambiguity is fixed by the additional condition (|9.60|) . Obviously solution H9.62|l obeys constraint (|9.60|l as well 
as the constraint H9.12|l . 



X. FUNCTIONAL Asig) IN THE CASE Nf ^2 
A. Functional AB{g,t) in the asymptotic regime t ^ 00 



In Sec. IV ill Dl we showed how the saddle point method can be used for the calculation of the functional Wo(g) 
(|8.25(l . The same saddle point method also allows us to compute the pre-exponential factor Asig)- To this aim we 
return to Eq. (|8.42l) and perform the integration over R more carefully, using Eq. (|4.22l) . As a result, we find 



27r 



exp [-Nc(t{h)] 



J dZD^^T^ (i?o-i) exp S^N, S^-tZ + i Indet ^ dyg{y)ye 



iZy 



where 



p-l 9cS 

— — = 
Vdetgeff 



ffeff = / dyg{y)ye'^y . 
lo 



(10.1) 



(10.2) 



(10.3) 



Next we perform the saddle point integration over Z which yields 
%i9) ' = 



— exp [~N^a{ii)] D j^T^ (R^ ) 



9^ , . 



dyg{y)ye'^y 



-1/2 



exp <j Nc { -iZ + - Indet 



dyg{y)ye'^y 



(10.4) 



where Z is determined by the saddle point equation 
instead of Z . Introducing the compact notation 



Lk — 



It is convenient to use the variable X = —iZ (|8.47f) 
dyg{y)y'+'e'''y , (10.5) 
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we can rewrite Eq. H10.2() in the form 



Using notation H1U.5|I . we find 

d 



Rn 



rtr 
^0 



Vdet Lo 



-1 + -Tr (L^^Li) 



Applying the identity 



we derive from Eq. H10.7|l 



9^ 



dZ 



In det 



Tr [L^^Li) = zTr - [L^^Li) 



-Tr 



Inserting Eq. pu.7|l into the saddle point equation (|8.4t)|) . we find 

Tr(Loiii) =2. 

This saddle point equation implicitly defines the dependence of X on 17 

X = X{g) . 

Inserting this dependence X{g) into Eq. Hl().5|l we see that Lk also become functionals of g only 

Lk = Lk{g) ■ 

Using Eqs. PU3|| . I^UT^ . and we find from Eq. ((Tini) 



$^(^) M^oo 2-^^+1 V2a(^)c^^. _ , ^,„, 



■ exp [-A'^cO-(m)] £',/3T3 



X {Tr - (Lq iLi)'] } exp 

Now we insert expression H4.8|l for ctn^ ^^nd apply the property of Wigncr functions (jA3l 



Nc I -iZ + -IndetLo 



{Tr[L^ii2- (io'^i)']} '^'exp 



Nc 1 -i^ + - IndetLo 



(10.6) 
(10.7) 
(10.8) 

(10.9) 

(10.10) 
(10.11) 
(10.12) 



In the leading order of the 1/Nc expansion, a{^) is given by expression 
0{N^) correction I^E to the energy NcE. Then 



1/4 



(10.13) 



\ V det Lq 

(10.14) 

We must also take into account the 



X(5) + - IndetLo 



X {Tr [l^^L2 - (io '^1)'] } exp |iV, 

This leads to the old result (|8.28() for the functional Was{g, t) 

W,sig,t)^Wo{g)-Et, 



'^NrJ V V det Lq 

1 



(10.15) 



(10.16) 



Wo{g) = Xig) + -ln[detLoig)] 



(10.17) 
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but now we also know the functional ABig,t) 

AV'i9:t) = cbV2TTTd^^j^ [Qo{g)]Ao{g,t) , (10.18) 
where cb is a constant independent of g and 

Qo{g) = -#=Vt ' (10-19) 
i/det Lois) 



Ao{g,t)=e^pi-tAE)[TT[L,\g)L2ig)~[L,\g)Li{g)]'}} . (10.20) 

Note that the functional AB{g,t) depends on the type of the state B. In the asymptotic regime t ^ oo this 
dependence appears in Eq. 1)10. 18|l only via the quantum numbers and via the constant cb- The origin of this 

simple dependence is obvious: the asymptotic distribution amplitude is controlled by the single operator with the 
lowest anomalous dimension corresponding to given TJ3T3. The result H10.18() can be used only for baryons obeying 
condition | J3I < T H8.29|l . because the underlying expression (|8.30(l for the asymptotic distribution amplitude was 
valid only for this case. For the lowest 0{N^^) baryon excitations with J = T (|3.1|l . the condition | J3I < T is satisfied 
automatically. 



B. Comments on the noncommutativity of limits Nc 00 and t — > 00 

In the previous section we have computed the functional AB{g,t) in the asymptotic regime. This calculation 
was based on the analysis of the double limit Nc ^ 00 and t ^ 00. As was discussed in Sec. IVIII these two 
limits do not commute. Therefore our results H10.18|) - (|10.20l) must be taken with a certain care. The calculation 
of Asigji) made in Sec. IX Al relied on the results of Sec. IVIII CI which were based on the explicit expression 
for the asymptotic distribution amplitude. Therefore the work of Sec. IX Al corresponded to the case when the 
asymptotic limit i — > 00 precedes the large- iV^ limit: the asymptotic limit selects the contribution corresponding to 
the lowest anomalous dimension (depending on the quantum numbers JJ3TT3), and only after that we use the large- 
Nc asymptotic distribution amplitude corresponding to these chosen quantum numbers JJ3TT3. Strictly speaking, 
the calculation of Sec. IX Al should be interpreted as a calculation of the Nf = 2 analog of the functional Aa{g) 1)9. 9|l . 
which appears in the problem of the diagonalization of the anomalous dimensions at large Nc- 

Although the method standing behind the calculation in Sec. IX Al is based on the limit 

firsts -> cx), then iVc 00 , (10.21) 

the form of presentation used in Sec. IX Al corresponds to the limit 

first A^c ^ 00, then i ^ 00. (10.22) 

Indeed, the functional AB{g,t) is defined only in the limit Nc 00. Before investigating the large-t behavior of 
AB{g,t) we must take the large- A'c limit in order to define this functional. 

To summarize, the results of the previous section were derived for the regime (|10.21(l . The possibihty of the 
extrapolation of these results to the region 110.22|l requires an additional investigation. 

If we stay in the "safe" regime (|10.21|) corresponding to the problem of the diagonalization of the anomalous 
dimensions, then we can rewrite Eq. (|10.18|l in the form similar to the Nf — 1 equation H9.59|l : 

ArnjAg) = V2TTTd?3J3 [Qo{g)] Ao{g) , (10.23) 

where 

io(.9) = {TT{L-\g)L,ig)-[L^\g)L^ig)]'}y'^' . (10.24) 

The zero mode factor \/2T + IZJ^^ [(3o(5)] in Eq- H10.23|I plays the same role as the Nf = 1 zero mode factor 
[^rot{g)V^ in Eq. (|9.59|) . Equation (|10.23|) contains no analogs of the nonzero mode factors £,^-,{g) of Eq. (|9.59|) 
because the representation (|10.23|l was derived for the asymptotic wave function corresponding to the lowest anomalous 
dimension (for given TT3J3). 
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The analogy between Eq. (|10.23|l and the Nj 
Nf = 2 zero mode Qo{g) obeys the same equation 



1 equations of Sec. IIX CI can be continued, 
as the Nf = 1 zero mode i^rot(ff): 



SWojg) ^ 6Qoig) 
Sg Sg 



For example, the 



(10.25) 



The check of this equation is straightforward but tedious. One has to use expressions (|1(J.17|I . (I1U.19|I for the functionals 
Wo{g), Qo{g) via functionals Lk{g) and the identity H1U.10|) for Lk{g). 



C. Zero modes and the structure of the functional Asig) 

In Sec. IIIIDI we have suggested the factorized form (|3.8f) for the functional Asig)- However, this naive factorized 
form is modified by zero modes. The precise form of this modification depends on the involved symmetries. For 
example, in the problem of the asymptotic distribution amplitude for Nf = 1 flavor we have only one nontrivial zero 
mode ^rot(5) l|9.49|l . This zero mode corresponds to the Abelian symmetry associated with the helicityJs. Therefore 
the zero mode factor [Crot(5)]''^ appearing in Eq. H9.59|l does not violate the naive factorized form (|3.8|l of the functional 
AB{g). 

The case of Nf — 2 quark flavors is more complicated. The results obtained in this paper give only indirect 
information about the functional Ab^q)- We know the functional As^g) 

1) in the toy quark model where it is given by expression (|4.12|) for the baryons with T = J, 

2) in the asymptotic i ^ oo case described by Eq. H1U.23|) . which was derived for baryons with the helicity J3 
constrained by the condition | J3I < T. 

Comparing Eqs. H4.12|l and H10.23() . we see that in both cases the functional Asisi) the form 

At=j,t,jA9) = const \Q{g)\ A^{g) . (10.26) 

In the toy quark model we have 

Q{g) ^ , ^0(5) = 1 (toy quark model) , (10.27) 

Vdet g 

whereas in the case of the asymptotic wave function {t 00) 

Q{g)^^M=, Aoig) = {TT{L^'ig)L2{g)-[L^'i9)Liig)]'}y'^' (asymptotic regime) . (10.28) 
VdetLo(.9) J J 

In both cases Q{g) is an SL{2, C) matrix. 

Comparing the results obtained in the naive quark model H1U.27|) and in the asymptotic limit H1U.28|) , we can make 
the conjecture that in the nonasymptotic regime the true QCD functional Axr^Jsig, t) for the lowest 0{N^^) excited 
baryons with J = T will have the same structure as Eq. 1)10.26(1 

AT=j,T,j,{9.t) = consti^i;^^^;^ [Qig.t)] Ao{g,t) , (10.29) 



det Q{g,t) = 1, 



(10.30) 



where Q{g, t) is some unknown functional of g with values in SL{2, C) which is determined by the dynamics of large- A'^j, 
QCD. 

Let us show that the conjecture 1(10. 29() and the SL{2, C) constraint 1(10. 30|l are compatible with the evolution 
equation. Applying the general evolution equation 16.20|l to the expression (|10.29|l . we see that it takes the form 



dP^rj/' [Qig,t)] 
dt D^-i[Q{g,t)] 



{g®g)-K- 



m{g,t) S Drjlf [Qi9,t)] 
Sg ^ Sg D^=J[Q{g,t)] 



(10.31) 



Let us show that this equation will hold automatically for any T — J,Tz,J^ and T' = J',T^,J^ if the functional 
Q{g,t) obeys the evolution equation 



d 

-g-^Qfs{g,t) = -\{g<E)g)-K^ 



SW{g,t) ^ 5Qfs{g,t) 



Sg 



5g 



(10.32) 
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Indeed, Eq. H10.32|l has an obvious property: If Q{g,t) obeys this equation, then any function F{Q) depending on 
the matrix elements Q wiU obey the same equation: 



Applying this property to the function 



we derive equation (|l().31|l . If we take 



then we find from Eq. (|10.33|l 



{g(g)g)-K' 



SW{g,t) ^ SF[Q{g,t)] 



5g 



5g 



F(Q) =detQ, 



d_ 

di 



detQ(5,t) 



{g ®g)-K- 



5W{g,t) 6detQ{g,t) 



Sg 



(10.33) 



(10.34) 



(10.35) 



(10.36) 



Obviously this equation is consistent with the SL{2, C) constraint l|l(J.30|l . 

Certainly the consistency of the conjecture H10.29|l . p0.30|) with the evolution equation and with the asymptotic 
case cannot replace the proof of this conjecture. In Rcf. 24| we present additional arguments in favor of Eq. (|10.29|) 
based on the large- iVc spin-flavor symmetry JJ], 2X, .22, _23J and on the soft-pion theorem for the baryon distribution 
amplitude 49] . In Ref. we also derive the correct version of the naive factorized ansatz (|3.8|) for the 0{N^) excited 
states, taking into account the modification of Eq. 1)3. 8() due to the zero modes. 



XI. CONCLUSIONS 



This paper shows how the systematic 1 /Nc expansion can be applied to the baryon wave function. The main idea of 
the method is to introduce a generating functional for the baryon wave function and to construct the 1 /Nc expansion 
for this generating functional rather than for the baryon wave function itself. 

A specific feature of this functional is its exponential dependence on Nc- The exponential part of this functional is 
universal for all low- lying baryons [including the 0{N~^) and 0{N^) excited resonances and meson-baryon scattering 
states]. The pre-exponential factors depend on the baryon (baryon-meson) state but these factors have a simple 
representation in terms of functionals associated with elementary excitations (or single meson scattering states). 

The developed formalism was applied to the light-cone baryon wave function (distribution amplitude) . It is shown 
that the structure of the 1/Nc expansion is compatible with the evolution equation. A nonlinear evolution equation 
is derived for the leading (exponential) term of the generating functional for the baryon distribution amplitude. The 
nonlinearity of this equation can be understood in terms of the well-known analogy between the large- A^c limit and 
the semiclassical approximation. The usual linear evolution equation is an analog of the quantum nonstationary 
Schrodinger equation, whereas the nonlinear evolution equation appearing in the large-A'c limit corresponds to the 
classical Hamilton- Jacobi equation. 

The nonlinear evolution equation derived in this paper was successfully used for the diagonalization of the anoma- 
lous dimensions of baryon operators of the leading twist. This problem was reduced to the linearized analysis of 
perturbations of the asymptotic solution of the Hamilton- Jacobi equation controlling the evolution in the large Nc 
limit. Although the nonlinear evolution equation has a rather nontrivial functional form, this problem was solved 
analytically, and simple expressions were found for the anomalous dimensions of baryon operators in two orders of 
the 1/Nc expansion. 
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APPENDIX A: BARYON WAVE FUNCTION IN THE TOY QUARK MODEL 

1. Wigner D functions 

The quark wave function 14. 2|) is represented as an integral over the SU{2) group containing Wigner functions 
^mim2 (■^) '^ith the weU known properties 

^^,^,(^1^2)- E Dl^rniRl)DL^,iR2), (Al) 



dR[Di^^„,^{R)] Dl,^,{R)^^-^S,,,S^,^,S^,^,. (A2) 
For the transposed matrix i?''' we have 

DL.mAR'') = D'in.mAR) ■ (A3) 

The characters of the irreducible SU (2) representations 

j 



have the properties 



Xi/2(i?) = Tri?, (A6) 



dRx*AR)XsAR) ^ . (A7) 



2. Normalization integral 

Let us compute the constant ctn^ appearing in Eq. 14.211 . This constant ctn^ is fixed by the normalization condition 
(El: 

\cTNfJ dR'J dRDl^^iR-') [Dl^,{R'-')r [Tr(i?i?'+)]^^ = 1. (A8) 

Changing the integration variable R RR' and using the properties of the D functions HA1|) - (|A2|) , we can simplify 
this integral: 

J dR'J dRDl^^R-') [Dl^M'-')V [Tr(i?i?'+)]"^= 

^ (2t\i)^ I dRT^^^miR'') (Tri?)"'^ . (A9) 
According to Eqs. (jA4|l - (|A6|) we have 

J^DlmiR-') (Tri?)^= = Xt{R) [xi/2iR)f' ■ (AlO) 
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Now we derive from Eqs. IjASp - IIA10|) 



where 



CTN^ = rj , (All) 



/at^j = j dRx*{R) (TvRf^ = J dRx*{R) [xi/2(i?)] • (A12) 

This integral has a simple meaning. It counts how many times the j representation appears in the decomposition of 
the tensor product of Nc spin-i representations. Combining this group interpretation of the In^,] with the standard 
rules of the spin addition, we arrive at the recursion relation 

with the initial conditions 

/i,i = l, /-i,^. =0. (A14) 
It is easy to check that the solution of Eqs. (|A13|I . (|A14|) is 



Inserting this result into Eq. l|Alip . we find 



3. Function (^TTsJa (5) 

According to Eq. H4.6|) 

<PTT,J, (<?) = CTN^ J dRDl^^ [Tr {Rg'^)] . (A17) 

In the case when g is an SU{2) matrix we can compute this integral changing the variable R^ R {g^'^) ^ 

^tt.jM = CTN^ j dRDl^Ag'^R'') (Tri?)'^^ 



T 

CTN. J2 I dRDl^,ig''^)D^,^,^{R-'){TTRr^ . (A18) 



Using the characters (|A4|) - ljA6|l . we can express this integral in terms of Ij.n^ I|A12|I 



j dRD^,^^{R-'){TrRf' = ^^ jdRY. Dl^{R-')i^^-Rf' 



= I dRxAR) {TrRf^ = ^TT^^'^^ ' ^^^^^ 

Thus 

^tt.jM - Dl,^ig-f-^^^ . (A20) 
Using the above result for ctn^ (|A11|I and the property of Wigner functions ljA3|l . we find 

^TT,j, {9) = Vi^Aj. (5) ■ (A21) 
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This result is derived for SU{2) matrices g. By definition H4.6() . 't'TT^Jaig) is an analytical function of g. Therefore 
we can "analytically continue" the equality HA21() to arbitrary SL{2, C) matrices g [assuming the SL{2, C) version of 
the Wigner function D'^^j^{g)]. 

Now let us consider arbitrary matrices g. According to the definition of i'TTiJiig) H4.6|l we have 



VTT'- 



Therefore 



0TT3J3(5) = (detg)^"/^0TT3,/3 5(det5) ^'"^ 



Since g (det g) is an SL{2, C) matrix, we can use the result (|A21|I 

gidetg)-'/'] - ^^h^M.j, Udetg)"^/' 



PTT3J3 



Combining Eqs. ljA23|l . (|A24|I and using the expression I|A11|I . we find 

ctn^ l 



(A22) 
(A23) 

(A24) 

(A25) 



APPENDIX B: EVOLUTION KERNEL 



A detailed discussion of the evolution equation, kernels and anomalous dimensions for the baryon distribution 
amplitude can be found in the original paper and in later publications, e.g. Here we list only those 

relations which are used in this paper. 



1. Nc dependence 



The leading-order evolution equations for the quark distribution amplitudes of mesons 



d N'^ — 1 a (ri) 

M^*mcson(2;i,a;2) = Ki2-^ l^^^^^^{xi , X2) 



and baryons 



M^<aryon(2;i,---,a;ArJ 



2Nr TT 



l<i<j<N^ 



(Bl) 



(B2) 



contain the same "pair interactions" Kij. The operator Kij acts on the variables Xi, Xj. The TVc-dependent factors in 
Eqs. ljBl|) and ljB2|) originate from the contraction of the color structure of the one-gluon exchange 



(B3) 



Sp {tH'') = ^6"'' 

with the color singlet projectors corresponding to the hadron states. This gives for mesons 



2Nr, 



-Sr.', 



(B4) 



(B5) 



and for baryons 



5Z *Cic'i*C2 4 ^4 



c'^C3...Cm^ ~ 2_/Y ^CiC2C3...CjVj. 



(B6) 
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2. Properties of the evolution kernel 

The evolution equations ljBl|l and ljB2|l are written in a compact form in terms of operators Kij . The action of K12 
is described by the kernel K^^'^'^ {xi,X2;yi,y2) introduced in Eq. H().4|l . The spin structure of this kernel is 

K'''^xuX2;yi,y2)^S''''^K+ixi,X2;yi,y2) + S''''''K~ixi,X2;yi,y2). (B7) 
We assume that the momentum conserving delta function is included into the definition of the evolution kernel: 

K'''^ixi,X2;yi,y2) ~ S{xi + X2 - yi - y2) ■ (B8) 
One can introduce homogeneous polynomials IsEj l 

B fr T 1 - V n\{n + 2)\{-lf , ^ 

Rn+l(.X„X2) - 2^^,(^^i),(„_^),(^_fc^^),^1^2 , (B9) 
fc— 



3 /2 

which can be expressed in terms of Gegenbauer polynomials C„ : 



n+l ^ „ c 

7n ^^^Z^fc (n + l)(n + 2) ^ ^ 



A:=2 



We have 



Here 



k=2 



for any function /. 
In particular, 



(BIO) 



Rn+l{xuX2) = -^{X,+X2TCI/^ f^^^\ 

n+l \xi+ X2J 

The polynomials (a;i, 3:2) diagonalize the evolution kernel: 

/ dyi I dy2K''^^^{xi,X2;yi,y2)Rn+l{yi,y2)yiy2 = ^7n''^Rri+lixi,X2)xiX2 ■ (Bll 

Jo Jo ^ 
The anomalous dimensions 7^^"^ are 



1*00 POO 

J dyi J dy2K^{xi,X2\yi,y2)Rn+i{yi,y2)yiy2 = -jltRn+i{xi,X2)xiX2. (B13) 



^s^s, ^^s,s.^+^^s,,-s,^-^ (B14) 



7+ = l + 4^i (B15) 



n+l ^ 2 



Since the kernel /v ^(xi, X2; yi, 2/2) contains the delta function (|B8|) . we immediately conclude that 

/ rfyi / dy2K^{xi,x2;yi,y2)Rn+i{yi,y2)yiy2f{yi + y2) = -^lnRn+i{xi,x2)xiX2fixi+x2) (B17) 



OC /"OO 

t^2/i / dy2K^{xi,x2;yi,y2)Rn+iiyi,y2)yiy2Siyi + y2- 
Jo 

--f^Rn+lixi,X2)xiX2S{xi + X2 - \) , (B18) 
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poo I'OO -|^ 

Jo Jo ^ 



Taking n = in these equations, we obtain 



OO /'OO 



Ay\ / c/?/2-fi^''"*^(a;i, a:^2;2;i, 2/2)2/12/2^(2/1 +2/2 - A) = -J^'''^xiX2(5(xi +X2 - A) , 
^0 ^ 



(B19) 



(B20) 



CXD /'OO 



dyx \ dy2K'''%xi,X2;yi,y2)yiy2e-''^y'+y'^ = -6^''' xiX2e'<^'+^''^ 



The evolution kernel can be decomposed in Gegenbauer polynomials 

00 - 

^Y.f(^n'\^y-^)cl'\2x~i) 

'e{y-x) 9{x-y) 



n=0 



1 

\x - y\ 



5.- 



where 



{l-x)y {l-y)x 
(n + 2)(n+ 1) 



(a; ^ 2/) , 



is the normalization constant for Gegenbauer polynomials: 



dz{l - z^)CU^ (z) (z) = /l„^„„ . 



(B21) 



(B22) 



(B23) 



(B24) 



Special care is needed for the singularities appearing &t x = y in Eq. ljB22|l . The precise integral form of this 
equality is 



dx 



where ^(a;) is an arbitrary function. 



00 g „i 

^Y.ir^n\'^y-i) / dxCl/H2x-i)x{i-x^{x) 

n=0 

(x) - (f>{y) 



(^) + 



F - y\ 



-9{y-x) + - — -9{x-y) 
V 1-2/ 



'-My) 



(B25) 



APPENDIX C: FUNCTION ALS Tks{g) 

In this appendix we study the properties of functionals Tks{g)- These functionals are defined by Eq. H9.22|l . They 
play an important role in Sec. IIX CI where we solve Eq. H7.11(l . 



1. Variational derivatives 

Let us start from the variational derivative of the functional M^sig) H9.21|l : 



Sgsiy) 

Combining this result with Eq. 19.22() . we find 



Sgsiy) Mos{g) 



5gs{y) 



(Ci) 



(C2) 



The variational derivative of the identity H9.45|l with respect to 5gs (y) yields 

^ 5Tis' (g) 



= 0. 



We insert Eq. into Eq. ifUS)! : 



Then 



where 



5X{g) _ N{g) 



5gs{y) Mos{g) 



[y-Tis{g)]ye 



-X{g)v 



N{9)^\Y.{T2A.9)-[T,Ag)f] 



2. Identities for the action of the evolution kernel K 

a. Action of K on 5Wo/ Sg ® SX/ 5g 
Now we want to prove the following identity 



{g<S>g)-K 

According to Eqs. H8.10|l and (|9.21() we have 



5Wo{g) ^ 5X{g) 



Sg Sg 



= 0. 



Using Eqs. (|C5|) and (|C8|I . we find 



SWojg) ^ ye-^(g)^ 
5gM 2Mos{g) ■ 

SWoig) ^ 6X{g) 



{g(E)g)-K 



5g Sg 



yie 



'X{g)yi 



2MosA9) 



N{g) 



E 



MosA9)MosM 



[{gs,<E>gsA-K'''' { 



yiy2e 



MosAy) 



'X{g)(yi+y2) 



[y2 - Tis, {g)] 



In order to compute the action of the operator K we use Eq. (IbT7|) 



5lS2 I 



K 



yiy2e 



-X{g){yi+y2) 



[y2 - TisAy)]] 



-X(g){yi+y2) 



yme 



-X(g){yi+y2) 



'Tis2{9) + ^ (2/1 +2/2) + ^(2/2 - 2/1) 



-Tis2{g) + ^ (2/1 + 2/2) 



+ ^7r^^(2/2-2/i) 



Inserting this into Eq. ljC9|l . we obtain 

'5Wo{9) ^ 5X{g) 



{g®g)-K 

N{g) 



Sg Sg 



N{g) 



E Wo'^' {-TuAg) + 1 [TuAg) + TuAg)]] - l^r^ [Tis 

SlS2 

Y^iYo^^^-YA'n[TisAg)-TisAg)] 
S1S2 



Here 



^"1^2 _ ^S1S2 _ S2S1 _ S2S1 
10 II — '0 II 



Tis, (5) - Ti,.^ {g) 



is symmetric in si, S2 whereas 
is antisymmetric. Therefore 

E (^0^^^ - ir^) [Tis, (9) - T,,2 (5)] - . 
S1S2 

Thus the RHS of Eq. !^Ju\ vanishes. This proves the identity ifUT)! . 

b. Action of K on 5Wo/5g (g) SMos/Sg 

Taking k = in Eq. HC1|) and using Eq. IIC7p . we find 

'SWoig) ^ SMosig)' 



ig^g)-K^ 



5g Sg 



Sss2y2e-'"^'^y' 



Inserting Eq. ljC8|l . we obtain 



{g(g)g)-K 



6Wo{g) ^ SMosig) 



Sg Sg 



= T.^^i9s.^gs)-K^ 



yiy2e 



-X(g){vi+y2) 



Here according to Eq. (|8.1|l 



Therefore 



Thus 



ViV2e 



-X(g){vi+y2) 



yiy2e 



-X(g){yi+V2) 



(gs, ® gs) ■ K''' 



yiy2e 



-X{g)iyi+y2) 



= -Ss,sMqs,Mos 



{g®g)-K 
As a consequence, we obtain the identity 

{g(S>g)-K 



SWAg) ^ SMosig) 



Sg Sg 
SWoig) ^ S Mo+ig) 



Sg SgMo-ig) 



:Mosig)- 



= 0. 
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We insert this expression into Eq. ljC22|l 



<5ff Sg 



1 



E 

Sl 

1 

4 



4Afo.i(g)A/o.(5) 



4A/o,,(g)A/o.(5) 



(ffsi I [(2/12/2) e-^(s)(^^+^^)] 



.m=0 



5*1 Tl— 



Thus relation H9.26|l is proved. 



3. Coefficients a^k and B^l^"^ 



a. Properties of a„k 



(C31) 



CoefHcients a„fe are defined by Eq. I)C23(I which can be rewritten in the form 

' = j2^n,cr{^)■ 



1 - Z 



Using the orthogonality of Gegenbauer polynomials (|B24|) . we find 



1 



ttnfe = T- / dz{l - z^) . „ 



1 - Z 



(C32) 



(C33) 



h. Properties of B'^^^^ 

Coefficients -6^4*^ ^^'^ defined by Eq. HC27() . According to Eq. (|ljl4() the anomalous dimensions 7^^*^ depend only 
on relative helicities: 



The same holds for the coefficients BH^'^ l[C27jl so that we can introduce B,.^ : 

^k7i - ^kn ' ^kn — ^kn ' 

In Sec. IIX CI we need the values for Bf.^ and BIq, 



Blk - - 



25, 



(fc + l)(fc + 2) k{k + l) 



[k > 0) 



(C34) 
(C35) 

(C36) 



Bio = 1 - 



k + 2 



fe+i , 
2E- (^^0)^ 



i=2 



which will be computed below. At k — we have 



-^00 - 



In order to derive relations (|C36|I and (|C37(I we rewrite Eq. (|C27p in the form 

.3/2 f 2/1-2/2 



n=0 



(C37) 



(C38) 



(C39) 
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Inserting Eq. ljC33|) . we find 



X / dz{l-z^) 



m=0 
^3/2, 



(C40) 



The summation over n on the LHS can be extended to infinity since Cn {z) is orthogonal to polynomials of degree 
smaller than n: 



(2/1 + y2f I dz 



1- z 



yi + y2 



k 



E BlmyTy2-"' ■ 



(C41) 



Changing the integration variable 



z^2x-l, 



(C42) 



we obtain 



^{yi + y2f / dxx{i-x) 
Jo 



fe+1 



H^K \yi + y2j 



ryiyi 



(C43) 



Taking (p{x) = (1 — x)'' in identity l|B25p and inserting the result into Eq. 1C43|I . we obtain 

\x-y\ 



(2/i+2/2r (l-y) -2 dx 





-9{y - x) + ^^e{x - y) 

y i-y 



" yi +y2 



(C44) 



Setting 2/1=0 here, we find 



BIq = 1-2 I dx 



2 



c . ^fc (1-2:) -1 

(1 - + ^ — - — 



k+1 



k + 2 



i=2 



(C45) 



Similarly taking ?/2 = 0, we obtain 

f-i 



-2 / da; 
/o 



<J..- (1 - xf + (1 - x)'^-^ 



26e 



(fc + l)(fc + 2) fc(fc + l)' 



(C46) 



Thus relations (IC36II and (IC37II are derived 
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